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ABSTRACT 


The Fleet Commanders-in-Chief often request recommendations from Commander, 
Mine Warfare Command, on schedules for naval mine transshipment. This thesis de- 
velops and implements a model, Scheduler for Mine Transshipment (SUMIT), that 
generates optimal schedules for transporting mines of a single type to suitable staging 
sites and for laying mines in mine fields. The model considers the number of available 
air, land and sea assets such as military aircraft, trucks, submarines and ships in finding 
optimal schedules for mine transshipment. SUMIT is designed to solve problems for 
scenarios in a limited region of the world that last several days and is based on inter- 
connected, time-expanded mine and mode networks. SUMIT is written in the General 
Algebraic Modeling System (GAMS) and is a mixed integer linear program in which all 
integer variables are binary. Ten realistic test problems are solved to demonstrate the 
viability of SUMIT and to compare the relative efficiencies of two model variants. One 


variant is on average 87% faster than the other. 


THESIS DISCLAIMER 


The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within 
the time available, to ensure that the programs are free of computational and logic er- 
rors, they cannot be considered validated. Any application of these programs without 
additional verification is at the risk of the user. 
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I. INTRODUCTION 


A. PURPOSE 

The Fleet Commanders-in-Chief often request recommendations from Commander, 
Mine Warfare Command (CMWC), who supports mining operations executed by the 
U.S. Navy, on schedules for naval mine transshipment in the event of war or for mine- 
laying exercises. Schedules for transporting mines to suitable staging sites and for laying 
mines in mine fields must be generated to meet the needs of mine warfare, given the 
number of available air, land, and sea assets (such as military aircraft, trucks, subma- 
rines, and ships). Only a limited number of military assets are capable of carrying mines 
due to special equipment requirements. Furthermore, since most mine transshipment 
and delivery assets are not solely dedicated to mine warfare, a mine transshipment 
scheduler must take into account both transshipment and delivery times and mine field 
and resource priorities. The current methods of planning mine transshipment schedules 
are manual, with no utilization of computerized optimization algorithms. Recognizing 
the need to computerize data and algorithms for many mine warfare problems, CM WC 
has initiated the installation of a computer system, the Geo-Operational Planning and 
Assessment System (GOPAS), to automate and optimize several aspects of mining op- 
erations [Ref. 1]. 

The purpose of this thesis is to build a computer-based mine transshipment model 
that can accept data stored by GOPAS and 1s automatic with respect to formulation and 
execution. While this thesis focuses on the rapid deployment of mines in regional real- 
time or exercise scenarios that take two to three days to complete, recommendations are 
also made for approaching problems of longer duration, such as a few months. The 
solutions of the Scheduler for Mine Transshipment (SUMIT) model can be used to re- 
commend not only feasible schedules for mine transshipment but, in manv cases, optimal 
(or close to optimal) schedules. While SUMIT recognizes time as a critical factor, its 
consideration of mine field and resource priorities allows the user to account for intan- 
gible aspects of the scenario by including both time and priority in the measure of ef- 


fectiveness (MOE) used in the model. 


B. BACKGROUND 
The mine transshipment problem is modeled, using two networks. Nodes in the 


network represent supply, transshipment, and demand sites. In the network for mine 


flow, arcs represent possible transportation routes between nodes in which modes, such 
as trucks, aircraft, ships, and submarines, carry mines. The network for mode flow is 
similar to the mine network, except that additional return arcs are needed to represent 
the return trips of modes back to their home base. 

]. Node Description 

Nodes can be separated into three broad categories: supply, transshipment and 
demand. The U.S. Navy has 13 Mobile Mine Assembly Groups (MOMAGs) located 
throughout the world which store nunes and function as supply nodes. These 
MOMAGs are typically collocated with large U.S. Navy and U.S. Air Force bases which 
have resources, such as aircraft capable of transporting and/or deploying mines. The 
mines at the MOMAQG Ss are maintained at six levels of readiness, from disassembled to 
fully-assembled (readv-to-deploy). If mines are disassembled, MOMAGs have the fa- 
cilities to prepare mines for delivery. The rate of assembling mines in preparation for 
deployment is called "build rate" and is measured in mines per hour. Some MOMAGs 
have the capability of setting up dual parallel assembly lines, thereby doubling the build 
rate. The U.S. Navy has an inventory of approximately 16 types of underwater mines 
and the build rate 1s known for each mine type [Ref. 2]. 

To decrease delivery times, fully-assembled mines are also pre-staged at lo- 
cations other than the MOMAGs. They can be stored at bases that have mine delivery 
air modes, on aircraft carriers and surface ships to be deployed via ship-based air modes, 
or on submarines to be deploved by the submarine itself. Currently, surface ships do 
not have the capability to deploy mines but a new class of mine-laying ships, called 
High-Volume Mine Layers (HVMLs), are planned to fill this gap. 

Transshipment nodes are typically U.S. Navy and U.S. Air Force bases with 
aircraft capable of supporting mine deployment. In this thesis, a military base that (1) 
has no pre-staged fully-assembled mines in inventory and (2) is not collocated with a 
MOMAG, is categorized as a transshipment node. 

CMWC has developed mine field plans for a variety of scenarios that could oc- 
cur throughout the world. In fact, the execution of all mine field plans simultaneously 
would deplete the entire supply of mines. Because the potential demand is much greater 
than the supply, all mine field plans are prioritized according to the importance of their 
objective. Approximately 90 percent of the time, mine field plans require only one type 
of mine and, therefore, generate only single-commodity scheduling problems. 

Node data that can be obtained from GOPAS to support the proposed model 
SUMIT include the number of each mine type stored at every MOMAG at each level 


of readiness, the number of each mine type pre-staged at given bases, the number of each 
mine type demanded at every mine field, the build rate of every MOMAQG for each mine 
type, the capability for a double assembly line, the priority of the mine field, and the type 
of node (land or sea). 

2. Mode Description 

Modes can be separated into two broad categories: transportation modes and 
delivery modes. Modes that cannot deploy mines in a mine field are categorized as 
transportation modes. These modes include land modes, such as trucks, and air modes, 
such as U.S. Air Force C-141 and C-130 cargo aircraft. Delivery modes can deploy 
mines and include air modes, such as U.S. Navy and U.S. Marine Corps shore-based and 
ship-based aircraft and U.S. Air Force B-52 bombers, and sea modes, such as subma- 
rnes. Due to weight and size constraints, most air modes and all land modes can only 
carry one mine type per trip. Sea modes, which will include HVML Ss in the future, are 
able to carry more than one nune type aboard. Because modes have multiple capabili- 
ties, they may also be critical to the success of other missions and can be prioritized ac- 
cording to the scenario. For example, if Anti-Submarine Warfare (ASW) plays a vital 
role in the scenario, a higher priority should be placed on non-ASW modes to deploy the 
mines. 

Since the mine transshipment problem includes mobile supply sites (aircraft 
carriers, ships, and submarines), differences arise in how mobile supply sites are incor- 
porated into the network. For this thesis, resources normally classified as modes are 
treated as nodes if they function as supply nodes and can only store mines. For example, 
aircraft carners and surface ships, which store mines but are not capable of laying them, 
are Classified as fixed nodes in the model. The aircraft stationed aboard the vessel 
function as modes that lay mines. To maintain consistency in the structure of the net- 
work, submarines and HVMLs, which are able to deploy mines, are treated as fixed 
nodes that have dummy modes stationed aboard to deliver mines. In addition, because 
aircraft carriers and ships mav not be located at the on-station point for launch of air- 
craft at the beginning of the problem, SUMIT allows the transit of aircraft carriers and 
ships at the beginning and end of the problem. This exception is not required for sub- 
marines which travel directly to the mine fields. 

Mode data supported by GOPAS include category (transshipment or delivery); 
dimension (land, air or sea); speed (nautical miles per hour); capacity for each mine type 
(mines per unit mode); time to load mines (hours); time to unload or deliver mines 


(hours); total tame to refuel, change crews, and conduct routine maintenance or repairs 


to prepare modes for their next trip (hours); maximum range (nautical miles); the num- 
ber of each mode available at each supply and transshipment node; and mode priority. 
3. Arc Description 

Because the flow of mines between nodes in the network is very structured, tlie 
onlv arc data that require support from GOPAS is the distance between nodes (nautical 
miles). Arc capacity (mines per trip) and length (time periods) can be calculated from 
the node data. Based on thc data input by the user, SUMIT only forms arcs for mine 
networks that are directed from supply nodes to supply nodes, supply nodes to trans- 
shipment nodes, supply nodes to demand nodes, and transshipment nodes to demand 
nodes. SUMIT generates the same arcs for mode networks, except the return arcs are 
also included and are directed back to the home base. [n Chapter II, the section entitled 
"Network Generation Rules" explains in detail other rules used to generate the net- 


Works. 


C. SCOPE 

Other researchers have studied the multi-commodity transshipment problem for 
mine warfare and for other military applications. Wingate and Zakary [Ref. 3] proposed 
a continuous variable model for multi-commodity transshipment problems, that could 
be applied to mine transshipment. However, their model was too general and did not 
address characteristics unique to the mine transshipment problem, such as the existence 
of modes (1.e., aircraft carriers and submarines) that also function as supply nodes. 

Collier, Lally, and Puntenney [Refs. 4, 5 , 6] developed continuous variable models 
for military deployment problems from the U.S. Transportation Command 
(IRANSCOM ) using sea and air assets. The mine transsliipment problem is smaller 
with potential supply sites limited to the 13 MOMAGs plus the bases and sea assets at 
which mines are pre-staged. The proposed model SUMIT is designed for regional 
problems spanning a time window of two to three days in which the total number of 
nodes can be about ten, while the TRANSCOM model developed by Puntenney covers 
movement requests among as many as 22 ports, planning general schedules lasting up 
to three months. Due to the importance of time in a short-duration mine transshipment 
problem, sea assets are limited to U.S. Navy vessels, located near the mine fields at the 
beginning of the problem, that either have the capability to deploy mines or have air 
assets aboard that are able to lay mines, while sea assets play a prominent role in ship- 
ping material in the TRANSCOM problem. Since the mine transshipment problem re- 


quires mines to be deployed as soon as possible, it does not include opportune delivery 


times, which were used to generate costs for the objective function in the TRANSCOM 
problem. Finally, data concerning limits on the mine loading and unloading capacities 
at non-MOMAG supply sites and restrictions on the number of aircraft allowed at 
transshipment sites at one time is not available in GOPAS. 

The scope of this thesis could encompass a multi-commodity transshipment model 
that optimizes schedules for deploying different mine types in a global scenario. Multi- 
commodity problems usually involve mixed-integer programs that use extensive amounts 
of computer resources to find optimal solutions for large networks. Because only one 
mine tvpe can be transported by a mode for each trip, the complexity of the mine 
transshipment problem increases. Since very little preliminary groundwork has been laid 
which specifically addresses schedule optimization for the mine transshipment problem, 
the scope of this thesis has been narrowed to a single-commodity approach. Because 
this thesis focuses on regional scenarios that require mine deployment to a couple of 
mine fields over a period of several days, the use of a single-commodity model is justified 
by the fact that (1) most mine fields require only one type of mine and (2) the majority 
of MOMAGs do not supply all 16 types of mines. Furthermore, the transshipment of 
mines over long distances may take more time than is allotted for the scenario, which 
then restricts mine supply to sites in the same region of the world as the mine fields. 
Thus, a global problem could conceivably be divided into several regional subproblems, 
which could be solved separately, if nodes contained in one regional subproblem were 


not contained in the other regional subproblems. 


D. OUTLINE 

This thesis is divided into four chapters. Chapter I is the introduction. Chapter II 
proposes and describes in detail two model formulations. Chapter III discusses the re- 
sults of ten test problems that compare the size and speed of the two model formulations 
and recommends procedures for executing SUMIT in the General Algebraic Modeling 
System (GAMS) [Ref. 7], a software package. Chapter IV contains conclusions about 
the methodology described in this thesis, discusses the weaknesses of SUMIT, and re- 


commends future enhancements of SUMIT. 


Ill. MODEL FORMULATION 


A. GENERAL DESCRIPTION 

Both proposed mine transshipment models contain two separate networks: one for 
the mines and another for the modes. The mine network is represented by a set of bal- 
ance equations that controls the flow of mines through all nodes. Two types of variables 
are contained in the balance equations for the mine network: one variable type repres- 
ents the shipment of mines from origin nodes to destination nodes via modes and the 
other represents the inventory of mines at nodes. The mode network is represented by 
another set of balance equations that controls the flow of modes to and from all nodes. 
The mine and mode networks are linked by a set of equations that relate the flow vari- 
ables in the two sets of balance equations. In the linking equations, sending mines via 
an arc contained in the mine network, forces flow on the corresponding arc in the mode 
network. In other words, mines cannot flow through the mine network unless there is 
a mode to carry them. Furthermore, the number of mines flowing through the network 
is limited by the capacities of the associated modes. 

The mine network is directed, meaning that arcs are ordered pairs of nodes. All 
paths must start at supply nodes, possibly flow through other supply nodes and trans- 
shipment nodes, and must arrive at demand nodes. Circuits, which are paths that start 
and end at the same node, are limited to transfers between supply nodes. Figure 1 il- 
lustrates a mine network, where s, 7, and d, respectively represent supply, transshipment, 
and demand at node i. Notice that the path from s, to 7, contains two arcs, which indi- 
cates that two different types of modes are available to transport the mines. [Ref. $] 

The structure of the mode network 1s an expansion of the mine network, which has 
additional arcs returning empty modes back to their origin node. The return arcs are 
critical because they prevent arcs that represent the same mode leaving at a different 
time period from being used before the mode has returned from its previous run. Figure 
2 depicts a mode network, associated with the mine network of Figure 2. 

Time-expansion of nodes in the two networks and in the linking equations is re- 
quired to allow modes to make multiple sequential trips in either transporting or laying 
mines. Different arc lengths (which are measured in time periods) also make time- 


expansion desirable. In essence, time-expansion expands every node over several periods 


Figure 1. Example of a Mine Network 


of time. The arcs of the time-expanded network connect earlier nodes to later nodes. 
Figure 3 illustrates a very simple time-expansion of a mine network. 

To take into account the build rate of disassembled mines, a MOMAG-supply node 
is split into two nodes separated by an arc. The MOMAG node (the origin node), which 
has no arcs entering it, assembles mines, and sends them to the supply node (the desti- 
nation node). Although dummy modes represent build rate to maintain consistency in 
the structure of the mine network, they are not included in the mode network since 
build-rate arcs can be used every period and there is no limit on availability over all time 
periods. The MOMAG node in the time-expanded mine network can have only one arc 
leave it per time period and this single arc must arrive at the supply node to which it is 
connected. Figure 4 illustrates a MOMAG-supply node split into two nodes: m, and 
SØ 

The following terms are defined to avoid misinterpretation of the problem formu- 


lation since they are not part of standard network terminology: 


TRANSPORTA TION 
LEG WITH LOADED 
MODE 


=>] 
RETURN LEG WITH 
UNLOADED MODE 





Figure 2. Example of a Mode Network 


e Mode group - a set of modes based at the same node that are grouped by their 
similarity or by their command structure, such as seven P-3 Orion aircraft based 
at Naval Air Station (NAS) Adak, Alaska, or a squadron of various aircraft capa- 
ble of carrying mines deployed aboard an aircraft carrier. 


* Mode unit - an element of the mode group, such as one P-3 Onion aircraft or a 
squadron (if the squadron is composed of different aircraft). 


e Run- around trip from a node to a destination via a mode group. 


e Mobile node - a node which is not fixed in position (such as aircraft carriers, sur- 
face ships and submarines). 


* On-station point - the position from which a naval vessel conducts operations. 


B. ASSUMPTIONS 
Several assumptions were made in formulating the model to find the optimal mine 
transshipment schedule. Most assumptions reduce the scope of the problem to decrease 


the size of the network. The model assumptions are: 





Figure 3. Example of a Time-Expanded Mine Network 


The model is single-commodity, i.e., only one mine type can be shipped. 


The mine network or the mode network need not be connected. For example, 
Atlantic-based mines will usually be deployed to mine fields in the Atlantic, 
Pacific-based mines will usually be laid in mine fields in the Pacific, and there need 
not be arcs (1.e., routes) connecting the mine and mode networks for the Atlantic 
fleet to the networks of the Pacific fleet. 


A mode group can have only one destination per run and cannot be split to start 
more than one run at different times. 


Multiple runs can be made by a mode group to different destinations. 


Two levels of mine readiness are assumed for mines stored at the MOMAG-supply 
nodes: disassembled and assembled. 


Mobile nodes that cannot deploy mines (such as aircraft carriers and surface ships) 
are treated as immobile nodes for the duration of the problem. The air modes 
stationed on the nodes deliver the mines. 


Mobile nodes that cannot deploy mines may transit to an on-station point at the 
beginning of the problem and to a new destination at the end of the problem. 
These transit distances must be input by the user. However, no additional transits 


DISASSEMBLED ASSEMBLED 





Figure 4. Example ofa MOMAG-Supply Node Split in a Mine Network 


may be made to reposition a mobile node at other on-station points. A mobile 
node can only lay mine fields at demand nodes accessible from the on-station point. 


e Mobile nodes that can deploy mines (such as submarines) are also treated as im- 
mobile nodes with dummy modes attached to them that function as delivery modes. 
The dummy modes maintain consistency in the structure of the model. 


* Mobile nodes cannot be resupplied. 
e Several mine fields located in the same area may be grouped into one demand node. 
e All modes must return to their origin node before the problem ends. 
* A mode that can either transport or deploy mines will always be classified as a de- 
livery mode to eliminate unnecessary arcs. 
C. FORMULATION 
l. Indices 
The first four of the five following indices are used in the formulations of both 
models, while the last index applies only to the second proposed model: 


e - an element of the set of all nodes in the orginal network, 


n 
N = (1,2,...,n,...,nx), where ny 1s the number of nodes. 


10 


e y - an element of the set of all modes in the original network 


NE a Where nais the number of modes. 

e p - an element of the sct of all time periods used for the time-expansion, 
PESO Dom vieres the number of time periods plus one. 

RN Cam clementior the set of integers used for iterative loops, | — (1, 2; OR 


where n, = nk. 
e g- an element of the set of incompatible arc groups, G = (1,2,...,8,..., Hc) , where 
Ng 1s the highest total number of groups, where 7, = np— 1 
The elements in the set P actually represent the start and end points of time periods 
considered by the model. lor example, the first period starts at p equals O and ends at 
p equals 1. The extra period starting at 7 — 1 and ending at np is needed to calculate 
constants for mine inventory after the last period of the problem. The procedure for 
determining incompatible arc groups 1s described in the section under “Generation of 
Incompatible Arc Groups.” 

The node and mode indices can be categorized by several subset indices repres- 
enting the function, dimension or arc position of the node or mode. For example, 
MOMAG nodes, which are always origin nodes, can only be positioned on land and 
transshipment modes travel on land or by air. This information is critical in constructing 
a network that has realistic arcs and meets the assumptions of the model. The following 
subindices are subsets of the node set N: 

e n1 - an element of the set of al MOMAG nodes MN, where me MNcN. 
e s-an element of the set of all supply nodes SN , where s e SNcN. 

e r-an element of the set of all transshipment nodes TN , where 1 e TNCN, 
e d-anelement of the set of all demand nodes DN , where de DNCN. 


s am element Ser the set of all «nodes “positioned on land LIN; where 
le LNE(MN U SN UTN). 


e c- an element of the set of all nodes located on or under the sea CN, where 
ce CNS(SN U DN). 


e ¡- an element of the set of all origin nodes IN for arcs in the original network, 
where te IN = (MN U SN U TN). 


e ;- an element of the set of all destination nodes JN for arcs in the original network, 
where je JN = (SN U TN Û DN). 
Figure 5 indicates relationships between the function subsets (NIN, SN, TN, and DN) 
and the dimension subsets (LN and CN). Figure 6 depicts the relationships between the 
function subsets and the arc position subsets (IN and JN). Notice that the unions of the 


four mutually-exclusive function subsets and of the dimension subset are identical to the 


node set N . Finally, the union of the arc position subsets also equals N but their 
intersection is not empty: (INNJN) = (SN U TN). The only node subset that the user 
must input in addition to the node set N is the sea node subset CN. SUMIT assumes 
all other nodes are land nodes. The proposed model determines the membership of the 
function subsets from the input data. 


The final group of subindices are the following subsets of the mode set M: 


e fi -an element of the set of all dummy build modes BM that transfer mines as- 
sembled at MOMAG nodes to adjacent supply nodes, where f e BMcM. 


e : - an element of the set of all transportation modes TM that can only transfer 
mines to transshipment or land supply nodes, where 7 e TMcMI. 


e ó- an element of the set of all delivery modes DM that can deploy mines into mine 
fields, where ó e DMCM. 


e À - an element of the set of all modes LM that travel on land, where 
À e LMS(BM U TM). 


* y - an element of the set of all dummy modes CM for sea nodes that can deplov 
mines into mine fields, where y e CME&DM . 


e c - an element of the set of all aircraft modes AM that transport or deliver mines, 
where a e AM<(TM LJ DM). 
The same dummy build mode can be used for all MOMAG nodes if they have the same 
build rate (the treatment of parallel assembly lines is mentioned in the next section). 
Figure 7 indicates the relationship between the function and dimension subsets of M. 
The function subsets (BM, TM, and DM) are mutually-exclusive sets whose union forms 
the mode set M. The union of the mutually-exclusive dimension sets (LM, CM, and 
AM) also is identical to the mode set M. In addition to listing the mode set M, the user 
must input all function subsets and must indicate both sea and land subsets. The pro- 
posed model assumes that all remaining modes are air modes. 
2. Data 
The following list of parameters is a brief description of data contained in both 
model formulations. All data used in both models is fully described in the section enti- 
tled "Detailed Data Description for Both Models." The first parameter 1s not actually 
contained in the computer implementation of the models but is created to simplify the 
presentation of the equations located in the sections of this chapter under "Model A" 


and "Model B." The value of first parameter, OBJ,,, is not directly input by the user 


ypu 


and is calculated by SUMIT from input data. 


e OBJ, the value contributed to the objective function by the mode flow for mode 


۱اا 


u leaving origin node / arriving at destination node j at time period p. 
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Figure 5. Relationship Between Node Function and Dimension Subsets 


e PAR,- the number of time periods with length TPER required to load mode y 
based at origin node i, travel to destination node j and unload (or deliver) mines. 


e PRT,, - the number of time periods with length TPER required to load mode ہر‎ 
based at origin node i|, travel to destination node j , unload (or deliver) mines, re- 
turn, and make ground preparations for turnaround. 


e AMT, - the supply or demand of mines at node 7. 


e XUP, -the upper bound on mine flow variables from origin node | to destination 


Un 


node j via mode p. 
e ZUP, - the upper bound on mines in inventory at node n. 


e MXT., -the maximum total time (h) that mode group „u can be absent from its or- 
igin node 7 in making all of its runs. 


e MXR, - the maximum number of runs that mode group u based at its origin node 
i is allowed to take. 


e RHX, - the value of the right-hand side of the mine flow balance equation for 
mines leaving node ^ at the beginning of time period p, where p is either the time 
period before inventory can start changing or the time period after inventory can 
stop changing. 
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Figure 6. Relationship Between Node Function and Arc Position Subsets 


e RHY,, - the value of the right-hand side of the mode flow balance equation for 


ipu 


mode pu P at origin node z at the beginning of time period p. 


mu 


Of the preceding parameters, only AMT,, MYT,,, and MAR, aze input by the user. The 
remaining parameters are derived from the input data. 
3. Decision Variables 
Five decision variables are included in the first proposed model, while only the 


first four decision variables are included in the second proposed model: 


* v- the objective value of the model, which is a measure of effectiveness that in- 
volves node and mode priorities and return time for arcs selected by the model. 


° - the number of mines sent from origin node i via mode 4 and arriving at des- 


X pu 
tination node j at the end of time period p, where x,,, 2 0. 


e J, - a binary variable that equals 1 Lo 0 and equals 0 otherwise. 


Upu 


e x,,- the number of mines remaining in inventory at node z at the end of time period 
Pp. 


Jp, - à binary variable that equals 1 if all of mode units in mode group yu remained 
at origin node / during time period p and 0 otherwise. 
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Figure 7. Relationship Between Mode Function and Dimension Subsets 


Mine flow and inventory is represented by positive, continuous variables. Since, for 
fixed values of y,,,, the resulting problem is a single-commodity network flow problem 
with integer supplies, demands, and bounds on flows, these continuous variables are 
guaranteed to be integer [Ref. 9]. The value of y,,, indicates whether or not mode p is 
used for the arc from i to j arriving at period p and does not reflect the number of mode 
units within mode group u that are actually needed to carry the mines. Likewise, Vou 
indicates whether or not the entire mode group remains at node i during time period p. 
In some cases, the model solution may imply that some, but not all mode units within 
a mode group, will be empty as they transit. Bv the assumptions of the model, these 
empty mode units cannot be diverted to meet other mine field demands. However, in 
reality, if the scenario warrants it (1.e., squadron integrity is not required for defense 
purposes), the empty mode units can remain at the origin node or be diverted to ac- 


complish missions, not related to the problem. 


4. Model A 
Model A, the first model proposed in this thesis, consists of the objective 
equation, flow balance equations, linking equations, other constraints and bounds on 
the variables, which are explained in the section after the formulation. XTE and YTE 
are two sets that represent the variables associated with the arcs that exist in the time- 


expanded networks. XTE is the set of all variables in the mine network, while YTE is 


—À À——À — 


Sree 


the set of all variables representing | flow from origin nodes to destination nodes in the 
mode network. Rules that explain membership in these two sets are explained in the 
section under "Flow Existence Rules." The time periods p! , p?, and p? are defined as 
follows for applicable equations in the model formulation: 

e p-—prcl. 

e pP SPT ARE 

e pp = p— PRT a, + PAR, 


mp 


Minimize 


Px NN 


- 2, T 2 2. OBI pp Spy (1) 


(| 0 23S, EYTE 


Subject to | شی‎ o fer 
Ny nx 
RR RO,  ر(‎ Y) Du REA, DNIE 
a ue nip'u inpu np 
j=l xə X pu eXTE i=] 4u9Xnpy ATE 
Ep! — = RHX, Vd ol (3) 
X dp! Xap i Xidpu KE dp P 


i=] H9 Xp, EXTE 


Ynp! ¥np T DE pu 7 فرظ‎ = REE وو‎ VM Po (4) 
| ! j^? نت‎ eYTE i» VER eYTE 
AUP jpu Viton Sinn 2 O0 IP O BINE (5) 
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PRI; 


DI SS MAT, Vin (6) 


i poy, “YTE 


> a J'üpu < MAR, X (7) 


=1 paYypou EVTE 


۷۶۰۰۰۰۰٠٠۱۱ ۶۷۷۷٣۳۹9 ۸۶۶٦٥ (8) 
PaE O VLJ Py (9) 
E = ZUP np (10) 
Viu € 10,1) Vi p,p (11) 


5. Comments about Model A 
The following list of comments explains the purpose of the equations given in 
the preceding model formulation: 
e Equation l is the equation for the objective value. 


e Equations 2 and 3 are mine flow balance equations for every node and time period 
in which mine flow 1s possible. Equation 2 1s for nondemand nodes and Equation 
3 for demand nodes. 


e Equation 4 is the mode flow balance equation for every origin node, mode and time 
period in which mode flow is possible. 


e Equation 5 links the mine and mode networks by forcing the mode flow variable 


Jy, tO lif the mine flow variable x,,,is positive. 


e [Equation 6 is optional and ensures that the total time that a given mode is away 
from a given base is less than MAT 


e Equation 7 is optional and ensures that the total number of runs for a given mode 
from a given base is less than NIIR. 


e Equations 8 - 11 ensure that all variables are bounded. 


In the computer code for the model, the equal signs are relaxed (1.e., are changed to in- 
equalities) in Equations 1 - 4 to make the model easier to solve. The equal sign in 
Equation | is changed to greater than or equal to (>), while the equal signs in Equations 


2 - 4 are changed to less than or equal to (x). 


6. Model B 

Model B is an alternate model proposed for SUMIT. Itis based on the premise 
that, by only allowing one mode arc in a group of incompatible mode arcs to be selected 
by the model, the model can monitor the movement of modes correctly and prevent a 
mode group from taking-off for its next run before the mode group has returned from 
its previous run, without using binary variables for mode inventory. The alternate model 
has one additional index (g) and more internal data, eliminates the need for the mode 
inventory variable y,,,, replaces Equation 4 with a single equation, and deletes Equation 
11. The alternate model makes the same assumptions as the original model. The net- 
work generation rules are also the same and the mine and mode flow existence rules are 
modified. The set, IG,, represents the incompatible arcs derived for group g. The gen- 
eration of the set IG, is described more fully in this chapter under "Generation of In- 


compatible Arc Groups." 
Minimize 
np 


ne Rs 
3 25 2 A 2 OBJ pu Jip (1) 


f=1j=1p=0 pay, €YTE 


Subject to 
HN : 
NES 2, Au > xw, 7 RHX, VneDN,p (2) 
D رھ‎ eXTE i=] 49 Xiu EXTE 
ee Y Do xg, RHX& Vd (3) 
i=l uD Xau EXTE 
Rx 
Vip = 1 Vi, H8 (12) 
J=1 p2 yy, EYTEMG, 
AUT T Jpn pu c Qo Vis rH Sex, E (5) 
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7. Comments about Model B 
All comments from Model A also apply to Model B, except for those that per- 
tain to Equations 4 and 11. The purpose of Equation 12, which replaces Equation 4 in 
Model A, is to ensure that only one arc in a group of incompatible arcs is selected. 
8. Detailed Data Description for Both Models 
The data discussed in the next five paragraphs is input by the user and must be 
manipulated by both models to generate the networks, calculate the data describing the 
networks, time-expand the nodes to form the time-expanded networks, and develop data 
for the objective and constraint equations. The input data in the following list is scalar. 
The first three scalars are required input and the last three are optional (default values 
will be assumed if no input is given): 


e TEND - the time (in hours, h) by which all modes transport or deliver mines, 
meeting all demand, and return to their origin nodes, where TEND > 0. 


e TPER - the length of the time period (h) used in the time expansion, where 
US IPER S TEND. 


e RELT - the value of the zero tolerance for comparing real data, where 
DES RELT ۱0ت‎ 007 Le WiIix—y[ —-REET,then wx = y 1s assumed. 


e MNLA - the minimum distance (in nautical miles, nm) that air modes can transport 
mines to prevent air modes from shipping mines over short routes intended for land 
modes, where MNLA z 0. The default is 0. 


e MXSS - the maximum distance (nm) allowed for supply to supply transfers, where 
MXSS > 0. The default is the maximum distance between nodes. 


19 


e DISC - a binary parameter in which 1 means to run the model without checking 
for disconnected networks and 0 means to check for disconnected networks and to 
run the model for the first connected network found. The default is 0. 


The model converts TEND into LTEND/TPER] time periods since it time-expands the 
network over time periods of length 7PER. If TPER 1s relatively small, the round-off 
error is less when converting times to time periods. However, a smaller time period 
length also yields more time periods for the problem, which expands the size of the 
model. Since a large model is harder to solve, a balance must be struck between TPER 
and the number of periods that will be formed. The user must also make TEND large 
enough to result in a feasible solution while making it small enough to cut down the size 
of the time-expansion. 

The data presented in the following list pertains to the node index n or any of 
its subindices. For all data pertaining to priorities in this thesis, a lower value implies a 
higher priority, e.g., a demand node with a priority of 1 1s more important than a node 
with priority 2. The purpose of the third parameter MD, is to allow the user to group 
mine fields that are close together into one node and to account for transit time needed 
to travel between mine fields. Since a mode group can only have one destination node 
per run, MXD, permits the mode group to travel to all mine fields within a demand node 
if it has the capacity to carry enough mines. The first parameter is required as input and 


the last five are optional: 


e AMT, - the supply or demand at node zr (mines), where, by convention, assembled 
supply and MOMAG disassembled supply amounts are positive and demand 
amounts are negative. The default is 0. 


e PRN, - the priority of node », where PRN, >0O. The default is the maximum node 
priority (or 1 if none are input). 


e MXD, - the maximum distance (nm) between mine fields within a demand node d, 
where MAD, > O. The default is 0. 


e DGD. - the distance (nm) that must be transited by the sea node c to reach its on- 
station point before launching aircraft to lay nunes, where BGD, z 0. The default 
value is 0. 


e EDD, - the distance (nm) that must be transited by the sea node c to reach a new 
destination by the last period of the problem after it completes its last run, where 
EDD, = 0. The default value is 0. 


e SPX,- the transit speed (nm/h) of the sea node, where SPX, > 0. If BGD, > 0 
or EDD, » 0, SPX, 1s no longer optional. 


The optional parameters BGD., EDD., and SPX, are not intended to be used for sub- 


marines and HVMLs since these modes transit directly to the mine field. 
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The following hst of input parameters describes mode characteristics indexed 
by mode or any of its subindices. The first seven parameters are required and the last 
one 1s optional: 


e CAP, - the capacity of one mode unit of mode pz (mines), where CAP, >0. CAP, 
is required for all modes except dununy build modes, whose capacities are com- 
puted by the model. 


e. MAL, - the maximum distance (nm) that mode 4 can travel from its origin node 
and be able to return, where MAL, » 0. MAL, is required for all modes except the 
dummy build modes. 


e SPD, - the average transit speed (nm'h) of mode uy , Where SPD, > 0. SPD, is re- 
quired for all modes, where the speed of a dummy build mode / is its build rate. 


e TLD, - the average amount of time (h) needed to load a mode group and take-off, 
where ED. 0. (he defaults 0. 


e TUL, - the average amount of time (h) needed to unload a mode group for trans- 
portation modes or deploy the mines for delivery mode groups, where TUL, > 0. 
The default is 0. 


e TGD, - the average amount of time (h) needed to spend on the ground after a run 
before reloading and taking off, where TGD, > O. The default is 0. 


e PRM, - the priority of mode jj, where PRM, 50. The default is the maximum 
mode priority (or 1 if none are input). 


The parameter CAP, can be used in two different ways. If a mode group is composed 
of the same type of mode unit, the capacity of a mode unit can be input as CAP,. But, 
1f the mode group is composed of variety of modes, CAP, should be the total number 
of mines that the mode group can carry. The user must take this distinction into ac- 
count when inputting the number of modes stationed at node. For example, if the mode 
group at a give node is composed of seven P3 Orion aircraft and the user inputs CAP, 
as the capacity of a mode unit, the number of modes at that node should be seven. 
However, if the mode group based at a given node is a tactical air squadron and the user 
inputs CAP, as the total capacity, the number of modes at that node should be one. 
For dummy sea modes based at mobile sea nodes that can deploy mines, CAP, should 
equal the supply of the node in which the dummy sea mode is stationed. The three 
ume-related parameters, TLD,, TUL,, and TGD, , are assumed to equal 0 for dummy 
build modes. 

The next list of parameters involves the mode groups based at their origin 
nodes. The first parameter is required and the last four parameters are optional: 


e NUM, - the number of mode units 4 in the mode group based at origin node i, 
where NUM, > 0. NUM, 2 lforatleast one mode group containing mode units 


wat every origin node /. The default is 0. 


e BGT,,- the earliest time (h) in which mode group u is available for loading at its 
origin node i for its first run, where 0 < BGT, < TEND. The default is 0. 


e EDT, - the latest time (h) by which mode group u must return to its origin node i 
after completing its last run, where O < EDT, <s TEND. The default is TEND. 


e MXT,, - the maximum total time (h) that mode group 4 can be absent from its or- 
igin node i in making all of its runs, where O < MXT, < TEND. The default is 
TEND. 


e MXR,, - the maximum number of runs that mode group u based at its origin node 
i 1s allowed to take, where O < MXR, € [TEND| TPER]. The default is 
PRENDE 
For dummy build modes, NUM,, should equal 2 if the MOMAG node m has dual par- 
allel assembly lines and 1 otherwise. This convention will double the build rate for par- 
allel assembly lines. If the number of parallel assembly lines at MOMAGs are increased 
in the future, then the user can account for this growth by setting NUM,, equal to the 
number of parallel assembly lines. For dummy sea modes stationed on mobile sea nodes 
that can deploy mines, NUM,, should equal 1. 
The following optional input parameter pertains to arcs in the mine and mode 
networks. It allows the user to rule out arcs in the networks that will be selected by the 
model: 


e XMP,, - a binary parameter in which 1 implies that no arc from origin node i| to 
destination node j via mode group 4 shall be allowed in the network. The default 
value is 0. 

Thus, by including XMP,, if necessary, the user can eliminate arcs that are not allowed 
because of obscure rules not accounted for in the model. 
The final input parameter is required and defines the distance between nodes: 


e DIS, -the distance (nm) between origin node i and destination node j. The default 
value is 0. 
Distances may be given in a table containing all distances between points. Distances 
may also be listed separately for arcs that will probably be used in the network. If the 
distance from node i to node j is input, then the model assumes that it also equals the 
distance from node j to node i. The former method of input ensures that the model will 
consider all possible node combinations for the networks, while the latter method is less 
tedious to enter if the user thinks only a few arcs are needed. 
Two important conventions must be followed to ensure the successful generation 
of the correct mine and mode networks. Entering "negative" distances allows the user to 


describe two properties of the networks without creating additional input parameters. 
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First, all distances between nodes that are on the same land mass must be less than 
—0.5 to indicate arcs that have potential land routes: DIS, < —0.5, where l # l'. Sec- 
ond, to identify a MOMAG node and its adjacent supply node, the distance between a 
MOMAG node and its adjacent supply node must be equal to —0.5: DIS,, = —0.5. 
This convention ensures that the MOMAG node is connected only to its adjacent supply 
node. | 

The next list of parameters is computed by the model, given the input data to 
describe the original and the time-expanded networks for mode flow and mine flow. 
These parameters assume that the network has already been generated. (Generation 
rules will be given in the next section of this chapter under "Network Generation 
Rules"). The model takes a conservative approach to rounding time periods up or down. 
For example, when converting time data to time periods, the model rounds up begin 
times input by the user (BGT,,) but rounds down end times (EDT,,). Rounding data in 
this way ensures that the time periods fall within time bounds input by the user. Like- 
wise, calculations of arrival time periods for one-way trips and return time periods for 


round trips are rounded up to be on the safe side. 


* NRLT - the value of the negative zero tolerance allowed when comparing real data, 
where NRLT = — RELT. 


e SUMA - the number of arcs in the mine network. 


èe ARC,- the set that indicates whether or not the arc from node i to node j via mode 


u belongs in the mine and modes networks. 
e ART), - a temporary set that is identical to ARC,,. 


e DNO,- a temporary parameter that serves two different purposes in the model: (1) 
to indicate when node vn has already been checked for calculating MXM, in one 
loop and for calculating MXS, in another loop and (2) to indicate that node n is 
included in the first connected network found. 


e NFC, -the set of nodes that does not form a component with at least one supply 
node and one demand node. 


o UA - the set that indicates whether or not the mine arc leaving node i and ar- 
riving at j at the end of time period p via mode p, belongs in the time-expanded 
network. 

e YTE,,, - the set that indicates whether or not the mode arc leaving node i| and ar- 
riving at j at the end of time period p via mode 4, belongs in the time-expanded 
network. 


e MNU,, - the set of all mine flow balance equations in which at least one arc leaves 
node n at the beginning of time period p or one arc arrives at node n at the end of 
time period p. 
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MDU, - the set of all mode flow balance equations in which at least one arc leaves 
node i at the beginning of time period p via mode yp or one arc arrives at node i at 


the end of time period p via mode u. 


RHX, - the value of the right-hand side of the mine flow balance equation for 
mines leaving node n at the beginning of time period p, where p is either the time 
period before inventory can start changing or the time period after inventory can 
stop changing and 


> ee min 0 2 ٦ if p is before. 
RHX,p 7 16 É 3 if pis after. 


RHY,, - the value of the right-hand side of the mode flow balance equation for 


ipu 


mode yw based at origin node n at the beginning of time period p , where RHY, 


tpu 


equals 1 if p is the time period before mode y can start its first run, equals —l if p 
is the time period after it can finish its last run, and equals 0 otherwise. Notice 
that, since mode flow and inventory variables are binary, 1 represents the presence 
of the supply of mode 4 at n when the problem starts and —1 represents the final 
return of the mode supply at the end of the problem. 


PAR), - the number of time periods with length TPER required to load mode yu 
based at origin node i, travel to destination node j and unload (or deliver) mines, 


A 
PAR = TPER 


if u e BM. 
+ (DIS, + MXD)) | (SPD, UP ET E iso 


PRT,, - the number of time periods with length TPER required to load mode u 
based at origin node i, travel to destination node j , unload (or deliver) mines, re- 
turn, and make ground preparations for turnaround, where 


TLD, + TUL, + TGD, 
PRTy, = TPER 


1 if u € BM. 
(2 |DISy| + MXD) | (SPD, TPER) otherwise. 


IND, - the indegree of (or the number of arcs directed into) node n. 
OTD, - the outdegree of (or the number of arcs directed out of) node nz. 


BGP,, - the earliest time period in which a run can start at origin node i and arrive 
at its destination node j via mode u, where BGP Z BOTT TIIR: 


EDP, - the latest time period in which a run can start at origin node i and arrive 
at its destination node j via mode û, where EDP,, € EDT,| TPER. 


iju 


MXM, - the maximum number of mines that could be sent through node n 1f all 
demand on paths containing node n were filled. 


MXS, - the maximum number of mines that could be sent through node n if all 
supply on paths containing node n were sent. 


OBJ,,,- the value contributed to the objective function by the mode flow for mode 
A leaving origin node i arriving at destination node Jj at time period p, where 


24 


ON LR IO TM CIDR 


Jijpu 
e ZUP, - the upper bound (mines) on mine inventory for node z, where 
ZUP, — min(MA M, , MXS,). 


e XUP,, - the upper bound on mine flow variables from origin node | to destination 


vn 


node j via mode u, where 

AUP,, = nuníZUP,, ZUP,, CAP, NUM, j. 
The parameters BGP,, and EDP,, are the keys to eliminating unnecessary arcs. For 
example, if the origin node i is a transshipment node, the model starts creating outgoing 
arcs after the earliest period in which mines could have arrived at node i. Likewise, the 
model does not create arcs that send mines to an origin node in periods after the last 
possible run could have left the origin node. This requirement cuts down on the number 
of continuous and binary variables, which increases the speed of model execution. 


The parameter OBJ 


gj, 1$ composed of two parts: the priority parameters and the 


time parameter. The priority parameters, PRN, PRN, and PRM,, enable the model to 
consider the relative priorities of nodes and modes during optimization. Because scaling 
problems may occur if the objective value grows too large, the user should only input 
relative priorities and the maximum priority should be less than five. Since the model is 
a minimization, the model tends to select mode flow across arcs for which the priority 
portion of the objective value is smaller (implying that the priority is higher). The sec- 
ond part of the objective value is the time parameter, p+ PRT,, . The model will again 
tend to select mode flow for which the time portion is smaller. Thus, mode flow that 
Starts sooner and has a faster round trip time will be considered more optimal. Because 
it 1s the product of the priority parameters and time parameters, the objective value of 
SUMIT incorporates both priority and time into a measure of effectiveness. 

Model B generates nine extra internal parameters and deletes one internal pa- 
rameter (MDU,,) from the original model. These additional parameters are needed to 
delineate the set of incompatible arcs IG, Additional parameters generated by Model 


B are: 


e IBEG -a scalar which is incremented by one in the loop that computes GBG,,, and 


GED,,, for arcs in the incompatible arc groups IG,, where 


IBEG > MNB,,. 


e MRN,- the minimum return period PRT, over all arcs from node | directed to 
node j via mode q. 


t2 
"1 


e MRX, - the maximum return period PRT,, over all arcs from node i directed to 
node j via mode 4. 


e MNA,, - the minimum arrival period PAR,, over all arcs from node i directed to 
node j via mode y for which PRT, = MNA, 


yu 
e MNB, -the minimum begin period BGP, over all arcs from node i directed to node 
J via mode q. 


e MXE,,- the maximum end period EDP,, over all arcs from node i directed to node 
j via mode p. 


e LSG, - the number of incompatibie arc groups leaving i via mode yp, where 
EDG’,, is EDG,, for which PAR, = MNA, and PRT,, = MRA, in the equation 


iju 
LSG, — min(EDG',) * 2- MNB,, — MRN, 


e GBG, - the group number g of the first incompatible arc grouping IG, in which 
the arc from i via mode u tọ j arriving at time period p appears. 


e GED,,, - the group number g of the last incompatible arc grouping IG, in which 


the arc from i| via mode yu to j arriving at time period p appears. 
9. Network Generation Rules for Both Models 
Before discussing the mine and mode flow existence sets XTE and YTE, rules 
for generating the original network are given since membership in the mine or mode 
networks is a criteria for membership in the flow existence sets. Initially, the set of eli- 
gible arcs is assumed to be the power set Nx NxM. Arcs that do not meet the rules 
for network generation are eliminated from the power set. The following list describes 
the existence requirements for arcs from origin node i to destination node Jj via mode بر‎ 


that belong to the original network: 


e An arc must leave origin node i and arrive at destination node j via mode u if 
NOM, > 0,0 4 j, DIS, #0; and jit SNP E 


e Any arc from i to Jj via p is eliminated if XMP; = 1. 


e Any arc between two supply nodes on land, ¿, je SNNILN , via a land transporta- 
tion mode, u e TM(|LM , is eliminated if DIS, > 0 or | D/S,| > MXSS. . 


e Any arc between two supply nodes on land, i,j € SN()LN , via an air transporta- 
tion mode 4 e TM(JAM is eliminated if D/JS,, « 0 and one of the following state- 
ments arê true: DIS |" = 1 SS" 


e Any arc between two supply nodes, ije SN, via mode yp is eliminated if 
1, JELN or ug TM. 


e Any arc is eliminated between two nodes i,j via land transportation mode 
مم‎ 218۸5۲۰٣۸۱1۰۷۸75 ٣ 


e Any arc is eliminated between two land nodes ije LN via air transportation mode 
He TMQAM if DIS, < 0 and | DIS, | > MNLA. 
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e Any arc from sea node c is eliminated if it does not arrive at demand node d. 


e Any arc directed to either supply node s or transshipment node z is eliminated if the 
mode is a delivery mode ó. 


e Any arc arriving at a demand node d via a nondelivery mode is eliminated. 


e Anv arc leaving a MOMAG node m and arriving at any node j for which 
DIS,; + —0.5 is elinunated. 


e Any arc entering a MOMAG node mt is eliminated. 


e Any arc leaving a transshipment node z and arriving at a supply node s is elimi- 
nated. 


e Any arc from /to/ via qu is eliminated rf its round trip time 7PER x PRT,, is 


yu 


greater than the total time allowed AfXT;, or the time difference between 
LODI and BGI: 


e Eliminate all arcs directed from node i if node í has no mines in inventory at the 
beginning of the problem and no arcs are directed into node i. 


e Any arc from i to j via uw is elininated if BGP, 2 EDP,,. 


e Any arc whose origin node / and destination node j is not contained in the frst 
connected component found. 


10. Generation of Incompatible Arc Groups 

The procedure for determining incompatible arc groups for Model B sets up in- 
compatible arc groups around an arc that has a return period equal to MRN, „and has 
an arrivai period equal to MNA, . Because such an arc has the minimum return period 
and the minimum arrival period for mode 4 leaving node |, the first run associated with 
this arc has the earliest arrival time period. The number of potential runs for this arc is 
also greater than (or equal to) the number of potential runs for other arcs that do not 
meet this criteria because arcs with greater return periods are forced to make fewer runs 
within the same amount of time. Figure 8 illustrates this point by depicting the re- 
lationship between arcs leaving the same node via the same mode arriving at different 
destinations with different arrival and return time periods. Focusing on this arc, then 
forms the greatest possible number of incompatible arc groups for mode p leaving i since 
it has the most runs and covers all possible time periods since its first run also has the 
earliest arrival time period. 

In the first step of the procedure, incompatible arcs groups are generated for 
arcs from i via mode p to j arriving at the end of time period p with return periods of 
MRN „ and with arrival periods of MRA,. Incompatible arcs for an arc from node i ar- 
riving at node Jj by the end of time period p are those that leave i before the mode group 


can return to i. Then, any arc that leaves node i after period p— MNA, and before 
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Figure 8. Relationship between Arcs from the Same Node to Different Destinations 


period p+ MRN,, — MNA,,, 18 not compatible. Adding MNA, to both sides implies that 
incompatible arcs arrive at node j after period p and before period p+ MRN, The 
procedure starts at IBEG = MNB., for the first group, iterates through a loop that cal- 
culates GBG,,, and GED 
the loop for group LSG,, is completed. GBG 


increments JBEG by one, and repeats the calculations unul 


pu ypu) 


„pu 1S then the minimum group number such 
that the following bounds are true for arrival time period p: 
IBEG < p < IBEG + MRN,. Likewise, GED,,, is the maximum group number for 
which the bounds are true. Figure 9 illustrates the formation of incompatible arc groups 
for arcs with minimum arrival and return periods using this procedure. 

The next step of this procedure then determines arc incompatibility for all re- 
maining arcs (which have return time periods greater than MAN, or arrival time periods 
greater than MNA,). It selects an incompatible arc by comparing the time period that 
the arc leaves node | and the time period in which it returns to node / with the take-off 
and return time periods of the incompatible arc groups formed in the first step. Once 
again, the procedure initializes IBEG to MNB, for the first group and iterates through 


the loop to calculate GBG,,, and GED, 


ypu ypu 


until group number LSG,, 1s checked. The 
bounds for period p are more complex in this step because the arrival time periods are 


not equal to MNA, or the return time periods are not equal to MRN, GBG, 1s the 


ijpu 
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Figure 9. Illustration of Step 1 of Generating Incompatible Arc Groups 


minimum group number such that the following bounds are true for arrival time period 


P: 


p < IBEG+(2x MRN,) — MNA,, + PRT;, — PAR 


yu 


Likewise, GED 


10 illustrates this step in generating incompatible arcs groups. 


y. 1$ the maximum group number for which the bounds hold true. Figure 

Although it appears that all incompatible groups have been generated, one more 
step 1s required because in some situations two arcs in IG, may be compatible. This 
situation occurs when the earliest arc in a group to leave for one destination 1s compat- 
ible with the latest arc in the same group to leave for another destination. Figure 11 il- 
lustrates the situation when two compatible arcs are in the same group. To avoid 
compatibility within incompatible arc groups, the compatible arc thar does not have the 


minimum return period PRT, 


yu And the minimum arrival period PAR,, 1s removed from that 
group. Removing this arc does not change its relationship with other incompatible arcs 
because it does not have the minimum return and arrival period and must have already 


been included in previous incompatible arc groups. 
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Figure 10. Illustration of Step 2 of Generating Incompatible Arc Groups 


The final step of the procedure is to eliminate the arcs added to IG, in the sec- 
ond step that are actually compatible with arcs placed in IG, in the first step. Looping 


through all incompatible arc groups IG,, where 1 < g < LSG, „ change GBG, p, tO the 


نام تح اللہ 


maximum of GBG,.,, and g + 1 if the following is true: 


inypu 
e PRT,, > MRN, or PRT,, » MRN, and 
e p> MNB, +(g- 1) + MRN, — MNA, + PAR 


yu" 
Change GED, „p, tO the minimum of GED, yp. and g — l if the following is true: 
6 PKI = MEN OLIM رب‎ and 


iju 


e p« MNB,* (g — 1) + MRN, — MNA, + PAR, — PRT, 


ur yu” 
Figure I1 illustrates the final step of the procedure that generates incompatible arc 
groups. 
11. Flow Existence Rules 

Once the mine and mode networks have been generated, membership in the 
time-expanded network for Model A is simple to express. Rules for membership in the 
existence sets are as follows: 
? Xp 


BGP 00ت‎ EDE 


yu 


e XTE if the arc from i to j via mode u is part of the original network and 
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Figure 11. Illustration of the Final Step in Generating Incompatible Arc Groups 


e y,, € YTE ifthe arc from i to j via mode yp is part of the original network and 
BGP,, S p SEDP,,. 


The mine and mode flow existence sets for the Model B are very similar to those 
of the original model: 


€ x,,€ XIE if the arc from i to j via mode yp is part of the original network and 


"Pu 


Den sp EDI 


e y,, € YTE if the arc from i to j via mode yw is part of the original network and 
BOP S PEED 


e y, € 1G, if y,,, e YTE and GBG,,, < g < GED 


Upu ypu ^ 
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HI. COMPUTATIONAL RESULTS 


A. COMPUTER IMPLEMENTATION OF PROPOSED MODELS 

GAMS, General Algebraic Modeling System, [Ref. 7] was selected for several rea- 
sons to implement the proposed models into computer code and to compare the relative 
efficiencies of the models. GAMS is a software package that incorporates 
FORTRAN -based solvers for optimization problems. The advantages of GAMS in de- 
veloping optimization models are that GAMS allows changes to be made simply, takes 
care of mundane details such as array sizes, generates all equations needed to solve the 
problem based on the algebraic expressions given in the model formulation, and uses 
relational databases to organize the data. Also, GAMS outputs the size of the problem 
in terms of the number of individual constraints (called "equations" by GAMS) and the 
number of discrete and continuous variables as well as information on the efficiency of 
the model in terms of the number of iterations and the computer resources needed to 
solve the problem. Such output can be used to choose the fastest model of the two 
proposed models for SUMIT. 


B. COMPARISON OF MODELS 

Since Model B eliminates the need for tracking mode inventory, it should have a 
little over half the number of binary variables generated by Model A in most cases. 
Fewer binary variables should make the alternate model easier to solve and, thus, faster 
than the original model. To test this hypothesis, ten different problems were solved by 
both models. Table 1 summarizes the input data for the ten problems. Since having 
more than three transshipment nodes for a regional problem is rare, the number of 
transshipment nodes is set at two for all ten problems. Also, the number of node-mode 
combinations, which is the number of positive NUM, parameters input by the user, is 
included rather than the number of modes. The number of node-mode combinations is 
more representative of the size of the model since all mode types are rarely stationed at 
the same origin node. To keep this thesis unclassified, the input data for the ten prob- 
lems do not represent actual data, but are realistic approximations of actual data. The 
resulting output concerning model efficiency is expressed in Table 2. Both models ar- 
rived at the same optimal solution for all ten problems, which demonstrates the validity 


of Model B. The problems were executed on an Amdahl 5990 mainframe using 


GAMS/ZOOM, where is ZOOM is one of the packages available with GAMS for mixed 


integer programs. 


Table I. INFORMATION ABOUT SCENARIOS FOR TEN TEST PROBLEMS 


Model End 
Number | time 













Number ofl Number of| Number of! Number of Number of 
Periods odes MOMAG/| Demand Node-Mode 
Supply Combinations 






po 
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Several important conclusions about the relative efficiencies of the two models can 
be made, based on the information presented in Table 2. For all ten problems, Model 
B produced fewer constraints than Model A, reducing the number of equations by an 
average of 15%. Both models have the same number of continuous variables because 
both models contain the same flow balance equations to track mine flow and mine in- 
ventory. Table 2 shows that, on average, Model B has fewer binary variables (by 39975), 
fewer iterations (by 87%), and smaller work parameters (by 76%). The work parameter 
controls the number of nodes waiting to be checked during the “branch and bound” 
portion of execution. If the work parameter is set too small, ZOOM will ternunate ex- 
ecution and ask the user to reset the work parameter. Thus, for problems similar to the 
test problems, the sixth column of Table 2 suggests that setting the work parameter at 
10,000 would suffice. In short, the data in Table 2 indicates that, for problems similar 
to the ten problems tested, Model B yields fewer constraints, binary variables, and iter- 


ations and smaller work parameters. 
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Table 2. RESULTS OF TEN TEST RUNS TO COMPARE BOTH MODELS 
Problem Number of | Number of | Number of | Number of Work 
and Model | Constraints | Continuous Discrete Iterations Parameter 


Numbers Variables Variables 
AE 54 | ae Sø 
A RE M — 413 | 100 
A 113 
A A .د ا‎ e. 
4A TS TT | 23,663 ٦ 
A os 
s | aa ——92| —  90| 06120] 3283. 
[ s |  n]| 6 TT 523 
عم‎ | — 1x7| — 96 — 84. 7608 
| éB |  —12| 96! 5|  3ens5| 233 
ا اع‎ ar ٥۲۵ قلله‎ 
| 6| s] %] 4“ 
Ds | | o sf sind) sa 
a o] ass 
os Pl xl s] 0 
]وی‎ 139] — 1| — 8| 187465 | SS 
Average BÎ IE — 1| — 2250| 12K 


Percent 15% 0% 39% 87% 16% 
Reduction 


Of the ten problems executed by Model A, problem nine had the highest number 





of iterations and required 12 minutes of CPU time to reach optimality. When problem 
nine was solved by Model B, CPU time was cut to 0.75 minutes. Problem seven was the 
most dramatic illustration of the reduction in CPU time, while problem three was the 


least dramatic. For problem seven, the CPU time for Model A was 3.5 minutes, which 
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was reduced to only 0.1 minutes when solved by Model B. The CPU time taken by 
Model A and by Model B to solve problem three was 7.5 minutes and two minutes re- 
spectively. The fact that Model B produces fewer constraints and fewer binary variables 
means that it requires fewer iterations and a smaller work parameter to reach optimality. 
This reduction in size and increase in speed, cut the CPU time by two-thirds in the 


Worst-case situation of the ten probleins tested. 


C. GAMS IMPLEMENTATION OF SUMIT 

The implementation of SUMIT and a sample problem in GAMS is given in Ap- 
pendix A. The output from the first run of the problem in Appendix A is contained in 
Appendix B. The input for and the output from the second, and final run, are contained 
in Appendices C and D respectively. Appendix E discusses several aspects of the GAMS 
implementation of SUMIT, including a description of SUMIT's output. 

Because large mixed integer problems are hard to solve, the user should limit the 
number of time periods to ten 1f the problem has more than ten nodes. However, since 
it 1s difficult to judge the end time and time period length for a new problem, the user 
should follow the recommendations given in Appendix E to find appropriate values for 
the end time and time period length. Notice that the efficiency data in Table 2 does not 
include the initial runs discussed in Appendix E because these runs are very quick. 

Other mixed integer program solvers are also available with GAMS. The solver, 
XA, solved the first run of sample problem in Appendix A in only 1$,859 iterations, 
while, as shown in Appendix B, ZOOM required 76,589 iterations. Thus, XA seems to 
be superior to ZOOM 1n solving mine transshipment problems using the SUMIT model. 

A way to attack large problems is to allow SUMIT to divide the large problem into 
several subproblems. SUMIT does not require that the nine or the mode network be 
connected. In other words, a path need not exist between all nodes m each network. 
If the network can be divided into multiple connected components, then SUMIT can 
split the problem into separate subproblems, one for each component. For networks 
that are large enough to tax the computer's disk space, running SUMIT for each con- 
nected component may solve a problem that could not otherwise be solved. Bv setting 
the scalar D/SC equal to 0, SUMIT will run the optimization on the first connected 
component of the network that it finds and will notify the user of the nodes not con- 
tained in the first component. The user must then rerun SUMIT for the nodes that were 


not in the first connected component. 
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The sample program given in Appendix A is an example of a large problem that was 
divided into subproblems by SUMIT. Appendix B contains the first solution output by 
SUMIT. SUMIT tells the user to rerun the program for nodes not included in the first 
component found. If there are at least one supply node and one demand node in this 
list that have not already been included in components already solved, the user should 
rerun SUMIT for the nodes indicated. The user may rerun SUMIT without deleting 
excessive amounts of input data by setting (1) NUM, equal to 0 for every node and 
mode combination that was included in the mine network for the first run of SUMIT 
and (2) AMT, equal to O for all demand already filled by the first run... Appendix C 
contains the input portion of SUMIT for the second run and illustrates this simple 
change. Notice that, since the second problem was smaller, the end time TEND and 
time period length T PER could be decreased: 

In many cases, finding the proven optimal solution is very expensive in computer 
time and memory, so GAMS provides an option that allows the user to specify an ac- 
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ceptable "optimality gap.” During the branch and bound portion of execution, ZOOM 
finds an upper bound and a lower bound on the objective value for the optimal solution. 
ZOOM terminates execution and reports an "INTEGER SOLUTION" if it has found 
a feasible integer solution in which the absolute distance between these two bounds di- 
vided by the lower bound is less than the optimality gap, called OPTCR. The integer 
solution reported when ZOOMS halts execution for this reason is not a proven optimal 
solution. In fact, when running the ten test problems using Model B, the solver found 
the optimal solution, without proving optimality, by using the GAMS default value of 
OPTCR (0.1) in six cases. Optunal solutions were also found, but not proven, for three 
of the remaining four problems by using an OPTCR of 0.01. Appendix E contains rec- 
ommendations for setting an appropriate OPTCR. 


D. REDUCTION OF MINE FLOW AND INVENTORY IN BOTH MODELS 

As stated previously, both models have the same number of continuous variables 
since they both contain the same flow balance equations for mine flow. Creating "rea- 
sonable" mine flow and inventory variables decreases the number of continuous vari- 
ables needed in the model. In the computer implementation of SUMIT, applying the 
network generation and flow existence rules ensures that only "reasonable" variables for 
nune flow and inventory are used. In the first portion of the “Pre-Model 
Manipulation" section of the program (see Appendix A), the set parameter ARC,, indi- 


cates the arcs that make sense for mine flow in the mine network. In the part of this 
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section preceding the generation of incompatible arc groups, the set parameter XTE,,, 
indicates the times for which mune flow associated with the arcs in ARC,, can exist. 
Likewise, the set parameter MNU, shows the only time periods in which potential 
changes in mine inventory may occur. Tables 3 and 4 illustrate the reduction of mine 
flow and inventory, respectively, for the ten test problems. They compare the potential 
number of variables (nune flow and inventory) with the number of variables actually 
created by SUMIT. Using the power set Nx N x M inflates the number of potential 
variables because the number of node-mode combinations can be estimated from the 
input data. Thus, the value for N x M is replaced by the number of node-mode combi- 
nations used earlier in Table 1. This swap is indicated by N-M in Tables 3 and 4 and 


provides a more accurate estimation of the number of potential variables. 


Table 3. MINE FLOW REDUCTION 
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The reduction in the number of variables for mine flow and inventory also affects 
the number of discrete variables because mine flow and mode flow are linked by the 
linking equations. In fact, the number of discrete variables (which is the number of 
variables for mode flow) created by Model B always equals the number of mine flow 
variables. Thus, by only including reasonable mine flow variables, SUMIT considers 
only logical mode flow variables. Decreasing the number of mode flow variables reduces 


the number of discrete variables, which increases the speed of execution for SUMIT. 
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Table 4. MINE INVENTORY REDUCTION 


Number N x N-M of MNUS Reduction 
Dai fo | | 
uw 59» 
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IV, CONCLUSIONS 


This thesis developed a single-commodity mine transshipment model, called 
SUMIT, for the rapid deployment of mines in a region of the world over a time window 
of two to three days. It considered two possible versions of SUMIT for computer im- 
plementation into GAMS. For ten problems tested, Model A and Model B arrived at 
the same optimal solution, but differed in relative efficiencies. For the size of the 
problems tested, Model B appeared to be faster and 1s, therefore, the model of choice for 
the implementation of SUMIT into GAMS. 

Transshipment models involving integer programs tend to be difficult to solve, even 
for small problems. When limited to a region of the world spanning a few days, the 
mine transshipment problem is small enough to attempt integer programming. This 
thesis explored a viable approach for solving small mine transshipment problems using 
integer programs by replacing the binary mode networks, which was done in Model A, 
with sets of incompatible arc groups for every mode, which was done in Model B. The 
results of the ten test problems in this thesis demonstrate the superiority of using in- 
compatible arc groups in the size and speed of the model for problem comparable to the 
test problems. 

In addition to providing optimal or nearly optimal schedules for mine transship- 
ment, SUMIT can be used as one of many tools for making decisions concerning mine 
warfare. SUMIT can give insight to logistical planners at CMWC in deciding the 
quantity and location of pre-staged nunes, the quantity of mines that should be fully 
assembled at MOMAGs, and the MOMAGs that should receive resources for dual as- 
sembly lines to decrease the amount of time needed to deploy sets of mine fields for the 
most likely scenarios. In preparation for mine deployment, SUMIT can be run for dif- 
ferent combinations of nodes and modes to indicate where modes should be based to 
decrease deployment time. Finally, SUMIT can be used to demonstrate the advantages 
of building new mine laying platforms, such as the HVMLs, bv not only illustrating the 
time that could be saved in mine deployment for short-term regional scenarios, but also 
indicating Which non-mine warfare platforms would be allowed to perform other critical 
missions. 

Several enhancements could be made to SUMIT to improve its user friendliness and 


expand the size of the problems that it can manage: 
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e Develop a "front end" program in another language to read user input from a data 
file and printout the GAMS program for execution. The front end program would 
read data from GOPAS and allow users to use SUMIT without directly interfacing 
with GAMS. 


e Develop a program in another language to read the GAMS output file and con- 
struct schedules that are more readable. 


e Allow a mode group from a given node to be split between destinations when ap- 
propriate. 


e Determine whether or not allowing more than two levels of mine readiness would 
be practical. This could be accomplished by splitting the MOMAG nodes so that 
each level of readiness for disassembled munes 1s represented by a MOMAG node 
and by putting one arc between MOMAG nodes for each level of readiness for 
disassembled mines and the associated supply node for assembled mines. 


[nclude constraints that take into account the limits placed on mine loading and‏ ٭ 
unloading capacities at non-MOMAQG supply sites and the restrictions on the‏ 
number of aircraft allowed at transshipment sites at one time. Currently, this data‏ 
is not available in GOPAS.‏ 


Finallv, the problem discussed in this thesis could be expanded for follow-on re- 


search in two ways: 


e Develop a model that encompasses world-wide scenarios of longer durations. This 
problem could be approached by adapting the integer rounding technique devel- 
oped by Puntenney to SUMIT [Ref. 6]. 


Develop a multi-commodity transshipment model bv applying the work of Colher,‏ ٭ 
Lally, and Puntenney [Refs. 4, 5 , 6].‏ 
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APPENDIX A. SAMPLE PROGRAMI, FIRST RUN 


STITLE * * * Scheduler for Mine Transshipment (SUMIT) * * * 
*----DOCUMENTATION-----------------------------------Ic-.-------.-.------- 


* Software to Optimize Schedules for Mine Transshipment 
* Developed for Commander, Mine Warfare Command 

* Written by Tammy L. Glaser, LT, USN 

* Complete documentation contained in thesis 


* 06 September 1991 

Eee | SANDS DOP BARSSCON FROFFOPT FONS FS FEE SEE SS ES SS SES SS SS SS DS S.e 
* Do not change the following options. 

SOFFUFPER OFFSYMXREF OFFSYMLIST 

OPTIONS SOLPRINT = OFF, LIMROW = 0, LIMCOL = 0; 

* Increase the following options only when recommended in the 

* solution report of the output. 

OPTIONS RESLIM = 10000, ITERLIM = 10000000, WORK = 10000; 

* Set OPTCR as recommended on page 72 of the thesis documentation. 
OPTIONS OPTCR = 0. 00001; 


MEAN TON SAND DADAS + Se == O a SS a و رت‎ ee ےا رت ا تا‎ 
* See Chapter 2 of thesis documentation to explain data in detail. 
SCALARS 
* The following scalars must be set to values greater than 0. 
TEND end time of problem (h) /48/ 
TPER length of time period (h) /6/ 
RELT zero tolerance for comparisons of real data / 0:0017 
* The following scalars are optional (enter O for default value). 
MNLA minimum distance for air modes over land routes (nm) 29517 


MXSS maximum distance for transfers between supply nodes (nm) /100/ 
DISC 1 to run entire problem and O to split into subproblems  /0/; 


* The first two sets must have the element O in them. Note that 

* N must contain all elements in CN(N), found in the optional section. 

* M must contain all elements in DM(M) (and TM(M), BM(M), CM(M), and 

* LM(M), which are found in the optional section). See page 7 of 

* thesis documentation to explain the MOMAG-supply node split. 

N nodes /0, M8GQ, S8GQ, MYRP, SORP; MIOJA, SIOJA; 

M11SC, S11SC, ANDERSNAFB, NASADAK, 
MCASIWKNI, NASJAX, CVNA, CVNB, SSN688A, 
SONGOBBB MINEPIELUIG HINEPIEED2 
MINEPIBDDS. MINEFIELDS/ 


M modes 7/0, TRUCK TCI O, 031212 P3 8525 TAGAIR, 
SSN688DM, DBRM/ 
DM(M) delivery modes /P3, B52, TACAIR, SSN688DM/ 


* P must start with O and end with the value of the number of time 
* periods plus one. G and I must start with 1. G must end with 

* the number of time periods, while I must end with the square of 

* the number of nodes (nonzero elements in N). 


P time periods /0 * 9/ 
G group numbers [187 
I iteration numbers /1 * 400/ 


* The following sets are optional. Enter 0 between slashes for sets 
* not needed for the problem. If CN(N) is empty (i.e., initialized to 
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* 0), nodes are assumed to be land nodes. If there are no transshipment 
* modes or disassembled mines, enter O between slashes for TM(M) or 

* BM(M). All modes are assumed to be air modes if CM(M) and LM(M) 

* are initialized to O. 


CN(N) sea nodes /CVNA, CVNB, SSN688A, SSN688B/ 
TM(M) transportation modes / TRUCK, C130, C141/ 
BM(M) dummy build rate modes / DBRM/ 
CM(M) dummy sea modes / SSN688DM/ 
LM(M) land modes [TRUCK / 
* Do not change the following evolution set. 
E mode evolution /ILD, TUL, TGD 


* Do not change the following set aliases: 
ALIAS (N,N1,N2,N3,N4); 

ALIAS CMMI MES 

7٦‏ ھ 

ALTAS (1511); 


PARAMETERS 
* The following parameter is required for node N: 
AMT(N) supply (+) or demand (-) at node N (mines) 
/M8GQ 200 
5860 25 
MORE 9 
1104 270 
S10JA SiS) 
M11SC 100 
S11SC 40 
ANDERSNAFB 89 
CVNA 161 
CVNB 100 
SSN688A 50 
SSN688B 40 
MINEFIELD1 Zu 
MINEFIELD2 -249 
MINEETELDO ZI 
MINEFIELD4 -201/ 


* The following parameters are optional (leave blank line with no 
* slashes for default values): 


PRNCN) priority of node N 
/MINEFIELDI 2 
MINEFIELD2 1 
MINEFIELD3 3 
MINEFIELDA 1/ 
MXD(N) maximum distance between mine fields at node N (nm) 
IMINEFIELD2Z 507 
BGD(N) transit distance for sea node N at beginning (nm) 
/ CVNA 33 
CVNB 754 
EDD(N) transit distance for sea node N at end (nm) 
/ CVNB 157 
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transit speed for sea node N (nm per h) 
25 


25% 


/CVNA 
CVNB 


SPX(N) 


The following parameters are required for mode M: 


in mode group M (mines per mode) 


capacity of mode unit 


/ TRUCK 
C130 20 
C141 20 
ES 8 
52 22 
TACAIR 51 
SSN688DM 507 
maximum distance that mode M may travel (nm) 
/ TRUCK 150 
C130 3000 
C141 3000 
۲3 1500 
B52 3000 
TACAIR 500 
SSN688DM 200/ 
transit speed of mode M (nm per h) 
/ TRUCK 30 
C130 380 
C141 410 
PS 350 
B52 480 
TACAIR 450 
SSN688DH 15 
DBRM 97 


CAP(M) 


MXL(M) 


SPD(M) 


The following parameter is optional (leave blank line for default): 


N CO a HA wma 


/ 


priority of mode M 
/ TRUCK 

C130 

C141 

ps 

B52 

TACAIR 

SSN688DM 

DBRM 


PRM(M) 


The following parameter is required for mode M at node N: 
number of mode M at node N (modes) 


PH WRAKF A DM a 


/M8GQ. DBRM 
S8GQ. TRUCK 
M9RP. DBRM 
S9RP. C130 
S9RP. P3 
M10JA. DBRM 
$10JA. P3 
M11SC. DBRM 
S11SC. TRUCK 
S11SC. C130 
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NUM(N,M) 


Y 


* 


لگ یی ہر ی بے بر یریم 


/ 


ANDERSNAFB. B52 
ANDERSNAFB. C130 
NASADAR. P3 
MCASIWKNI. TACAIR 
NAS JAX. P3 

CVNA. TACAIR 
CVNB. TACAIR 
SSN688A. SSN688DM 
SSN688B. SSN688DM 


* The following parameters are optional (leave blank line for default): 
earliest beginning time for mode M at node N (h) 
/SSN688A. SSN688DM 
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latest end time for mode M at node N (h) 
/ CVNA. TACAIR 


24/ 


maximum total time mode M can be absent from node N (h) 


maximum number of runs for mode M from node N 
/ ANDERSNAFB. B52 
ANDERSNAFB. C130 
SIRE C150 


2 
2 


2/ 


l means eliminate runs from nodes N to Nl via mode M 
/ S9RP. MINEFIELD2. P3 14 


S10JA 
DES 
1157 

138 
75 
uS 
0 
1028 
1028 
9 
2681 
1983 
71032 
849 
7293 
2144 
1388 
2170 
IOS 
15205 
7414 


M10JA 
1157 
1157 

738 
75 
0 
2055 
1028 
7028 
1172 
2881 
1083 
1032 
849 
7295 
2144 
7388 
2170 
1015 
1203 
7414 


BGT(N,M) 


EDT(N,M) 


MXT(N,M) 


MXR(N,M) 


XMP(N,N1,M) 


٠ 
3 


* The following two parameters are required. 
* READ PAGE 22 OF THESIS DOCUMENTATION, WHICH EXPLAINS WHEN NEGATIVE 
* DISTANCES MUST RD: 
TABLE DIS(N,N1) distance between node N and node N1 (nm) 


SORP 
1345 
1345 
05 
0 
738 
138 
7754 
7754 
1364 
3614 
1824 
7768 
1580 
8020 
259 
6125 
2999 
1725 
097 
8152 


M9RP 
1345 
1345 
0 
2053 
138 
738 
7754 
7754 
1364 
3614 
1821 
7768 
1580 
8020 
2581 
6125 
2909 
1725 
897 
8152 


S8GQ 
20.5 
0 
1345 
1345 
1157 
151557 
7501 
7501 
250 
2:439 
1651 
7485 
1566 
7789 
2405 
7812 
2459 
1819 
2207 
nono 


M8GQ 
0 


- 5 
1345 
1345 
1157 
1157 
7501 
7501 

=30 
2939 
1631 
7485 
1566 
7789 
2405 
78 12 
2459 
1819 
2207 
7797 


111580 511806 ANDERSNAFB NASADAK MCASIWKNI 


Do SM 


2999 


290 


44 


7301 


501 


M8GQ 

S86Q 

M9RP 

S9RP 

M10JA 
S10JA 
M11SC 
S1156 
ANDERSNAFB 
NASADAK 
MCASIWKNI 
NASJAX 
CVNA 

CVNB 
SSN688A 
SSN688B 
MINEFIELD1 
۶ك‎ 
MINEFIELDS 
MINEFIELD4 


+ 
M8GQ 


2 9 1631 
3614 1821 
3614 1821 
2881 1083 
2881 1083 
4222 5942 
4222 5942 
2939 1653 
0 1824 
1824 0 
4209 0 
2088 267 
4511 5942 
T 1067 
4546 6305 
744 2053 
2032 952 
3850 2-053 
4551 6305 
SSN688A  SSN688B 
2405 7812 
2405 7812 
2881 8152 
2881 8152 
2144 7414 
2144 7414 
4899 420 
4899 420 
2421 703] 
777 4546 
1067 6305 
4900 366 
1324 6546 
5164 S77 
0 9 
5249 0 
110 5220 
1244 6413 
3073 8023 
25270 174 
MINEFIELDS MINEF1ELD4 
2207 7597 
2207 7797 
897 8152 
897 8152 
1203 7414 
1203 7414 
7607 548 
7607 548 
2188 7812 
3850 455] 
2053 6332 
7660 463 
1786 6578 


230 
1364 
1364 
1172 
1172 
7520 
7520 

0 
2939 
1653 
7500 
1563 
7801 
2421 
7831 
2455 
1800 
2188 
7812 


CVNB 
7709 
7789 
8020 
8020 
723 
4293 
294 
294 
7801 
7801 
6216 
35 
6445 
0 
5164 
227 
5146 
6294 
7817 
500 


110857102 
1519 
1819 
1725 
1725 
1015 
1015 
6032 
6032 
1800 
2052 

352 
6051 
234 


7501 
7754 
7754 
7028 
7028 
20:5 
0 
7520 
4222 
5942 
-124 
6175 
294 
4899 
420 
4866 
6032 
7607 
548 


CVNA 
1566 
1566 
1580 
1580 
849 
849 
6175 
6175 
1583 
2088 
267 
6188 
0 
6445 
1324 
6546 
1344 
234 
1786 
6578 


ELD1 
2459 
2459 
2909 
2909 
2170 
2170 
4866 
4866 
2455 

745 
1090 
4869 
1344 


7501 
7754 
7754 
7028 
7028 
0 
Oss 
7520 
4222 
5942 
-124 
6175 
294 
4899 
420 
4866 
6032 
7607 
548 


NASJAX 
7485 
7485 
7768 
7768 
7032 
7032 
-124 
-124 
7500 
4209 
5350 

0 
6188 
366 
4900 
366 
4869 
6051 
7660 
463 


MINE ST 


S86Q 

MORE 

S9RP 

M10JA 
S10JA 
M11SC 
S1156 
ANDERSNAFB 
NASADAK 
MCASIWKNI 
NASJAX 
CVNA 

CVNB 
SSN688A 
SSN688B 
MINEFIELD1 
MINEFIELD2 
HINEFIELD3 
MINEFIELD4 


4 
M8GQ 

S8GQ 

MORE 

S9RP 

M10JA 
S10JA 
11150 
10ت‎ 
ANDERSNAFB 
NASADAK 
MCASIWKNI 
NASJAX 
CVNA 

CVNB 
SSN688A 
SSN688B 
MINEFIELD1 
MINEFIELD2 
MINEFIELD3 
MINEFIELD4 


4 
M8GQ 

S8GQ 

MYRP 

S9RP 
M10JA 
S10JA 
M11SC 
6 6 
ANDERSNAFB 
NASADAK 
MCASIWKNI 
NASJAX 
CVNA 


CVNB 5146 6294 7817 500 


SSN688A 110 1244 3073 5270 
SSN688B 5220 6413 8023 174 
MINEFIELD! 0 1265 3081 5258 
7272٦ 1265 0 1817 6459 
MINEFIELDS 3081 1817 0 8109 
MINEFIELD4 5243 6459 8109 0; 
TABLE TEV(M,E) time needed to complete evolution E for mode M (h) 

TLD TUL TGD 
TRUCK 1. 00 2. 00 1. 00 
C130 4. 00 2.00 8. 00 
CRM 4. 00 2. 00 8. 00 
P3 1. 00 0. 25 4.00 
Boz 4.00 0.50 8. 00 
TACAIR 1. 00 0:25 2. 00 
SSN688DM 0. 00 1:50 0.00; 
eee AD A A A A A À A A A A A AA A A A A A A A A AA RK KK KKK 
پو‎ * 
* Input is complete. Do not change any statements after this line * 
* Ds 
yc cv eve vec ves ee deve eve eee yes eye dee ye eese deve ve eve ye eve ve veo seve He Fe ve se se ste ا‎ e dee dee desee ve e eee eve eee پا از‎ Kee 
*----PRE-MODEL MANIPULATIONS----------------------------.-..----.-.------ 
SCALAR 


NRLT negative zero tolerance for comparisons of real data 

STOP stopping criteria for loop that finds connected component 
SUMA number of arcs in the mine network 

IBEG counter that is incremented in incompatible arc group loop; 


SETS 
NO(N) dummy node /0/ 
MN(CN) momag nodes 
SN(N) supply nodes 
TN(N) transshipment nodes 
DN(N) demand nodes 
LN(N) land nodes 
MO(M) dummy mode /0/ 
AMCM) air modes 
ARG CNN Lat) arcs in network from node N to node Nl via mode M 
ART(N,N1,M) temporary version of ARC 
DNO(N) arcs in first connected component found 
NFC(N) nodes that do not form a complete component 
XTE(N,N1,P,M) arcs in time-expanded network for period P 
MNU(N,P) mine flow balance equation used for N at period P; 


* Verify correctness of the inputted data. 

ABORTS(RELT LE 0) "ERROR: RELT must be positive." , RELT; 

ABORTS(RELT GE 0.00100001) "ERROR: RELT must be less than 0.001.", RELT; 

ABORTS(TEND LE RELT) "ERROR: TEND must be positive.'', TEND; 

ABORTS(TPER LE RELT) "ERROR: TPER must be positive.'', TPER; 

ABORTS(TPER GE (TEND + RELT)) "ERROR: TPER must be less than TEND.", 
TEER: TEND; 

ABORTS((MNLA + RELT) LE 0) "ERROR: MNLA must be nonnegative. '', MNLA; 

ABORTS((MXSS + RELT) LE 0) "ERROR: MXSS must be nonnegative." , MXSS; 

ABORTS((ABS(DISC - 1) GE RELT) AND (ABS(DISC) GE RELT)) 
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"ERROR: DISC must be equal to O or 1.", DISC; 


* Initialize or set default values for certain parameters. 

TEND = FLOOR(TEND / TPER); 

NRLT = -1 * RELT; 

MXSSS(MXSS LE NRLT) = SMAX((N,N1), ABSCDIS(N,N1))) + (2 * RELT); 

AMT(N) = FLOOR(AMT(N))SCAMT(N) GE 0.5) + CEIL(AMT(N))SCAMT(N) LE -0.5); 
MXL(M)SBM(M) = -0. 5; 

NUM(N,M) = FLOOR(NUM(N,M)); 

DS NOI, ٭‎ NOC(N) * MO(M)) = 

DIS(N,N1)S(NO(N) * NO(N1)) = 6; 

DIS(N,NI)SCABS(DIS(N,N1)) LE RELT) = DISD, 


TM(M)SMO(M) = NO; 
BM(M)SMO(M) - NO; 
CN(N)SNO(N) = NO; 
CM(M)$MO(M) = NO; 
LM(M)SMO(M) = NO; 


* Verify the correctness of the inputted data. 
ABORTS(SUM(MSDM(M), 1) LE 0.5) "ERROR: No delivery modes entered”, DM; 
ABORTSCSUMCNSCAMTCN) GE 0.5); AMT(NJ) LE O. >) 

"ERROR: No supply amounts (+) entered.", AMT; 
ABORTS(SUMCNSCAMI(N) LE -0.5), ABS(AMT(N))) LE 0.5) 

"ERROR: No demand amounts (-) entered.", AMT; 
ABORTS(SUM(N, AMT(N)) LE -0.5) 

"ERROR: Total demand (-) exceeds total supply (+) AMT; 
ABORTS(SMIN(M$(NOT (MO(M) + BM(M))), MXL(M)) LE RELT) 

ERROR: All MXL(M) must be positive. 

"Did you forget the maximum range for dummy sea modes?'", MXL, CM; 
ABORTS(SMINCHS(NOT (MO(M) + BMOD)), SPD(M)) LE RELT) 

“ERROR: All SPD(M) must be positive. 

"Did you forget the maximum range for dummy sea modes?”, SPD, CM; 
ABORTSCSMINCCN,M), NUMC(N,M)) LE -0.5) 

"ERROR: NUM(N,M) must all be nonnegative.", NUM; 
ABOBTSCSUMCC(N,N1);, ABSCDISCN,NI))) LE RELT) 

"ERROR: No distances (DIS(N,N1)) inputted", DIS; 
ABORTS(SUM((N,N1,M)SCCABS(XMP(N,N1,M) - 1) GE RELT) AND 

(ABS(XMP(N,N1,M)) GE RELT)), 1) GE 0.5) 

"ERROR: All XMP(N,N1,M) must be equal to 0 or 1.", XMP; 


* Categorize nodes and modes into appropriate subsets. 
MNCNOS="YESSCCAMTOØN) (GE 0.5) AND CSUSCHSBNCHO S; NUMCN;HO) GE 0.5955 
SN(N) 2» YESS((NOT (NOCN) -* MN(N)))SC(CAMT(N) GE 0.5) 

OR (CCABSCAMT(N)) LE 0.5) AND (SUM(NISMN(NI1), 

۲۰٣۸۳5۲0150٦1076 0,5) LE REEI)) CE 0.5)))): 


TN(N) 2 YESS((NOT (NO(N) + SN(N)))SCABSCAMT(N)) LE 0.5)); 
DN(N) = YESS(AMT(N) LE -0.5); 

CN(N) = YESS(CN(N) + DN(N)); 

LN(N) = YESS(NOT (CN(N) + NO(N))); 

AM(M) = YESS(NOT (LM(M) + CM(M) + BM(M) + MO(M))); 


* Include all logical arcs in network. 

ARC(N,N1,M)SCCCMN(N) + SN(N) + TN(N)) * (SN(N1) + TNCN1I) + DN(N1))) 
SCCNUM(N,M) GE 0.5) AND (ABS(ORD(N) - ORD(N1)) GE 0.5) 
AND ((ABS(MXL(M)) - ABS(DIS(N,N1))) GE NRLT) 


AND (ABS(DIS(N,N1)) GE RELT))) = YES; 
* Eliminate all arcs as indicated by user. 
ARCC(N,NI,M)SCABSCXMPCN,N1, M9) IU LE TELT) = NO; 
* Eliminate all arcs between supply nodes on land via land modes 
* that have nonnegative distances or have distances greater than MXSS. 
ARC(N,N1,M)SCCSN(N) * SNCN1) * LN(N) * LN(N1) * TM(M) * LM(M)) 
S(CDIS(N,N1) GE NRLT) OR ((MXSS - ABS(DIS(N,N1))) LE RELT))) = NO; 
* Eliminate all arcs between supply nodes on land via air trans- 
* portation when the arcs are over land and absolute distances are 
* less than MNLA or their absolute distances are greater than MXSS. 
ARC(N,N1,M)SC(C(SN(N) * SNCN1) * LN(N) * LN(N1) * TM(M) * AM(M)) 
SCCCCMNLA - ABS(DIS(N,N1))) GE RELT) AND (DIS(N,N1) LE NRLT)) 
OR ((MXSS - ABS(DIS(N,N1))) LE PELT))) = NO; 
* Eliminate all arcs between two supply nodes in which at least 
* one node is on land or that use nontransportation modes. 
ARC(N,N1,M)S(SN(N) * SN(N1) * (NOT (LN(N) * LN(N1) * TM(M)))) = NO; 
* Eliminate all arcs between nodes via land transportation modes 
* for which the distance is nonnegative. 
ARC(N,N1,M)SCCTM(M) * LM(M))SCDISC(N,N1) GE NRLT)) = NO; 
* Eliminate all arcs between land nodes via air transportation 
* modes when the arcs are over land and absolute distances are 
* less than MNLA. 
ARC(N,N1,M)SCCLN(N) * LN(N1) * TM(M) * AM(M)) 
SCCCMNLA - ABS(DIS(N,N1))) GE RELT) AND (DIS(N,N1) LE NRLT))) = NO; 
* Eliminate all arcs from sea nodes to nondemand nodes. 
ARC(N,N1,M)SCCNCN) * (NOT DN(N1))) = NO; 
* Eliminate all arcs to either supply or transshipment nodes 
* via delivery modes. 
ARC(N,N1,M)SC(CSN(N1) * TN(N1)) * DM(M)) = NO; 
* Eliminate all arcs to demand nodes via nondelivery modes. 
ARC(N,N1,M)SCDNCN1) * (NOT DM(M))) = NO; 
* Eliminate all arcs leaving MOMAG nodes to nodes for which 
* the distance is not equal to -0.5. 
ARC(N,N1,M)SCMNCN)SCABS(DISCN,N1) + 0.5) GE RELT)) = NO; 
* Eliminate all arcs arriving at MOMAG nodes. 
ARC(N,N1,M)$MN(N1) = NO; 
* Eliminate all arcs from transshipment nodes to supply nodes. 
ARC(N,N1,M)SCTN(N) * SNC(N1)) = NO; 


* Verify that only one arc leaves each MOMAG node. 
ABORISCSMAXCNSMNCNO, SUMCCNAIGMOSARGCGONSNISGEDOSUIODAGESU SS) 
"ERROR: At least one MOMAG node has more than one arc leaving it. ", 
"MOMAG nodes and the current network is listed as follows: ", MN, ARC; 


PARAMETERS 
PAR(N,N1,M) time periods to travel from N to N1 via M one-way 
PRTENSN1SM) time periods to travel from N to Nl via M round trip 
IND(N) in degree of (number of arcs entering) node N 
OTD(N) out degree of (number of arcs leaving) node N 
BGP(N,N1,M) earliest time to include arc from N to Nl via M 
EDPCNSNT.M) latest time to include arc from N to N1 via M 


MXM(N) maximum number of mines through N to fill demands 
MXS(N) maximum number of mines through N from supply nodes 
ZUP(N) upper bound on mine inventory at N 

XUP(N,N1,M) upper bound on mine flow from N to N1 via M 
RHX(N,P) fixed value on mine inventory at N at period P 
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via M 
via M 


return period PRT from N 
return period PRT from N 


arrival period PAR from N via M if 0 


ia M 
M 


last group number needed from N via M 


begin period BGP from N v 
end period EDP from N via 


m 
m 
m 
m 
m 


MRN(N,M) minimu 
MRX(N,M) maximu 
MNACN,M) minimu 
MNB(N,M) minimu 
MXE(N,M) maximu 
LSG(N,M) 


GBGUN.NISP.HMo first eroup nümber from N to-NI via MN at P 
GED(N,N1,P,M) last group number from N to Nl via M at P; 


Initialize or set default values for certain parameters. 


X(N1, PRN(N1))); 
X(M1, PRM(M1))); 
R); 
TEND; 
TPER * TEND; 
ben SPOR ۰ 
STOP, P; 
to 'TEND'.", TEND, G; 
CON STOP 
STOP, I; 


*X 


* Convert times to time periods. 


SMA 
SMA 


MAX(1, 
MAX(1, 


0(N1) * MOCM)) = 
O(N1) * MO(M)) = 
oo aan 1); 
= FLOOR(MXT(N,M) / TPE 


N 
N 


CAP(M) = FLOOR(CAP(M)); 
PRN(N)SNO(N) = 
PRM(M)SMO(M) - 


PRN(N)S( (NOT NO(N) )SCPRN(N) LE RELT)) 
PRM(M)S((NOT MO(M))$(PRM(M) LE RELT)) 


MXD(N)S(NOT DN(N)) = 0; 
BGD(N)SNO(N) = 0; 
EDD(N)SNO(N) = 0; 
SPX(N)SNO(N) = 0; 
BGT(N,M)S(NO(N)*MO(M)) 
EDT(N,M)$(NO(N)*MO(M)) 
SUMA = 0; 
MXT(N,M)$(NO(N)*MO(M)) 
MXR(N,M)S(NO(N)*MO(M)) 
GBG(N,N1,P,M)S(NO(N) * 
GED(N,N1,P,M)S(NO(N) * 


1 


MXR(N,M)SCMXR(N,M) GE RELT) = 
MXT(N,M)SCMXT(N,M) GE RELT) 


MXT(N,M)SCC(NOT (NO(N) + MO(M))) 


SCCMATCN,N) LE RELT) AND (ONUMC(N,M) GE 0.5))) 


EDTCN,MOSCCNOT (NOCN) -F MO(M))) 


Calculate capacity for dummy build rate modes. 
= FLOOR(SPD(M) * TPER); 


Verify the correctness of the inputted data. 


S((EDT(N,M) LE RELT) AND (NUM(N,M) GE 0.5))) 


x 
CAP(M)SBM(M) 


Ye 


Slory = TEND. + 1, 


ABORTS((ABS(CARD(P) - 2 - TEND)) GE RELT) 


contain integers from 0 
the value orf SIE, 
GE RELT) 

contain integers from 1l 


GE RELI) 

contain integers from 1 
the value of STOP ).' 
O REG LEO) 

all be nonnegative. 
RELA EEO) 

all be nonnegative.  , MXD; 
TFERELT ) LE O) 

all be nonnegative. , BGD; 
ERELT) -LE 0) 


tt PRN: 
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The set P must 

"(see the next line for 
ABORTS(ABS(CARD(G) - TEND) 
The set G must 
(CARD(N) - 1) ** 2; 
ABORTS(ABS(CARD(I) - STOP) 
The set I must 
"(see the next line for 


"ERROR: 


"ERROR: 
STOP = 


"ERROR: 


ABORTS((SMIN(N, PRN(N)) 
"ERROR: | PRN(N) must 
ABORTS((SMIN(N, MXD(N)) 


"ERROR: | MXD(N) must 
ABORTS((SMIN(N, BGD(N)) 
"ERROR:  BGD(N) must 
ABORTS((SMIN(N, EDD(N)) 


"ERROR:  EDD(N) must all be nonnegative.", EDD; 
ABORTS(CSMIN(N, SPXCN)) + RELT) LERO) 

"ERROR: SPX(N) must all be nonnegative.", SPX; 
ABORTS(SUM(NSCCCBGD(N) GE RELT) OR (EDD(N) GE RELT)) 

AND (SPX(N) LE RELT)), ٦ 

"ERROR: If BGD(N) or EDD(N) is inputted, SPX(N) must be positive." 

BGD,- EDD; SPA: 
ABORTS(S! MIN(MS(NOT MO(CM)), CAPD TEU 2), 

“ERROR: CAP(M) must all be positive. 

"Did you forget the capacity for dummy sea modes?", CAP, CM; 
ABORTSCCSMINCM, PRHMCMOO s RELD)SDESOD 

"ERROR: PRM(M) must all be nonnegative.", PRM; 
ABORTS((SMIN((N,M), BGT(N,M)) + RELT) LE 0) 

"ERROR: BGT(N,M) must all be nonnegative.", BGT; 
ABORTSCCSHIN( CNM), EDICN MO) ۶۱ + 1 71 

"ERROR: EDT(N,M) must all be nonnegative.", EDT; 
ABORTS((SMIN((N,M), MXT(N,M)) + RELT) LE 0) 

"ERROR: MXT(N,M) must all be nonnegative.'', MXT; 
ABORTS((SMIN((N,M), MXR(N,M)) * RELT) LE 0) 

"ERROR:  MXR(N,M) must all be nonnegative.", MXR; 
ABORTS((SMIN(N, SPX(N)) LE 0) 

"ERROR: SPX(N) must all be nonnegative.", SPX; 
ABORTS((SMIN((M,E), TEV(M,E)) + RELT) LE 0) 

"ERROR: TEV(M,E) must all be nonnegative.", TEV; 
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* Adjust BGT and EDT if nonzero BGD and EDD is inputted. Convert 

* BGT and EDT to time periods. 

BGT(N,M)S(NOT (NO(N) + MO(M))) 
= CEIL(((BGD(N) / SPX(N))$((SPX(N) GE RELT) AND (NUM(N,M) GE 0.5)) 
٭‎ ۶ 

EDT(N,M)SCNOT (NOCN) + 4+ ۵۷| 
= FLOOR(((EDD(N) / SPX(N))$((SPX(N) GE RELT) AND (NUM(N,M) GE 0.5)) 
٭‎ AE 


* Calculate PAR taking into account load, travel, and unload/ 
* delivery times. Calculate PRT taking into account load, 
* travel, unload/delivery, return travel, and ground times. 
PAR(N,N1 ,H)SCARCCN,NI, M)SCSPD(M) GE RELT)) 
= CEIL(((TEV(M, TLD') + TEV(M,'TUL')) / TPER) + 1$BM(M) 
+ (ABSCDIS(N NIDDA MXD(N1)) / -CIPBER * 50) ۹ ٘“ 7۶۳ 
PRI(N,N1,M)SCARC(N, N1 4 ) ۷١ 
= CEIL(((TEV(M,'TLD') + TEV(M,'TUL') + TEV(M,'TGD')) / TPER) 
+ (((2 * ABS(DIS(N,N1))) + 0) ۰)۹ ۱ٔ e 
+۰۱۹) (٦ 


* Eliminate all arcs from N to Nl via M if PRT is greater than 
* MXT or (EDT = BGT). If node N has no supply and no incoming arcs, 
* eliminate all arcs directed from node N. 
ARC(N,N1,M)SCARC(N,N1,M)SCCCMXT(N,M) - PRT(N,N1,M)) LE NRLT) 

OR ((PRTCN,N1,M) - CEDT(N,M) - BGTCN;MOODOGE 0 - NG; 
ARC(N,N1,M)S((SUM((N2,M1)SARC(N2,N,M1), 1) LE 0.5) 

AND (ABS(AMT(N)) LE 0.5)) = NO; 
ART(N,N1,M) = ARC(N,N1,M); 


* In a loop, take the following steps until a loop in which 
* no arcs are deleted from the mine network is made: 
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* Calculate SUNA, IND; and OID. Initialize DNO ando MAM. 


* Calculate EDP and MKM by checking all possible paths 

* in the network. Reinitialize DNO and MXS. Calculate 

* BGP and MXS by checking all possible paths in the network. 
* (In these loops, DNO indicates whether or not MXM or MXS 
* has already been updated for origin node of the arc being 
* checked. ) Calculate EDP for arcs leaving MOMAG nodes. 

* Eliminate arcs from N to NI via M for which PRIT is 

* greater than (EDP - BGP). 

* If node N has no supply and no incoming arcs from 

* other nodes, eliminate all arcs directed from node N. 

* If DISC equals O, then set 

* DNO to 'yes' for all nodes in the first connected 

* component found. If DISC equals 1, then set DNO to 'yes' 
* for all nodes (except dummy node 0). Eliminate arcs whose 
* nodes are not contained in the first component found. 


LOOP(I1S(ABS(SUMA - SUM((N,N1,M)SART(N,N1,M), 1)) GE RELT), 
SUMA = SUM((N,N1,M)SART(N,N1,M), 1); 


IND(N) = SUM((NI1,M)SART(NI,N,M), 1); 

OTD(N) = SUM((N1,M)SART(N,N1,M), 1); 

DNO(N) = YES; 

MXM(N) = 0; 

LOOP(NS(DN(N)SCCIND(N) + OTD(N)) GE 0.5)), 
DNO(N) = NO; 


LOOPCCNI.MOSARICNIN, M) 
EDP(N1,N,M) = MINCEDT(N1,M), TEND) - PRT(N1,N,M) + PAR(N1,N,M); 
MXM(N1)S(DNO(N1) + DNO(N)) = ABSCAMT(N)) * MXM(N1); 
DNO(N1) = NO; 
LOOP( (N2 MIJSART(N2 NIMD, 
EDEBCNZ NIM) =OMAXCEDPONZ,NI,MI), 
HINCEDPCNISNSMOSS PARCNI,N,;MO, MINCEDICN2,H1), TEND) 
- PRT(N2,N1,M1) + PAR(N2,N1,M1))); 
MXMCN2)SCDNOCN2) + DNO(N)) = ABSCAMT(N)) + MXM(N2); 
DNO(N2) = NO; 
LOOPCCNS MZ SARICNS N2 02), 
BDECNS NZ 2 = MAAC CEDP(N3; N2, M2); 
MINCEDPENZSN MU) PARCN? NISMI, 
MINCEDICN3 M2), TEND) 
ے٣‎ ۴7٦٥۷3٦۱79120706 ۸۲)۷3 ۳21220 
MXM(N3)S(DNO(N3) + DNO(N)) = ABSCAMT(N)) + MXM(N3); 
DNO(N3) = NO; 
ZOOP(CN% M3)SART(N4,N3,M3), 
EDP (NG NS MOI S BANCEDPCNASNO MO. 
BINCEDPCNS;N2.H2) - PARCNS NO Mo, 
MINCEDT(N4,M3), TEND) 
PRI 10۰5 ]13( PARONG,NS MS 
MXM(N4)S(DNO(N4) + DNO(N)) = ABS(AMT(N)) + MXMCN4); 
DNO(N4) = NO)))); 
DNDONIASCNOT DNCNY)) 2 YES): 
BOPONSNI,MISARTON,NI SM) 
SEPARCN NESCII DIF LE 0.5) + TENDSCINDON) CE 0.5) 
SSBGTONSEOS 
DNO(N) = YES; 
MXS(N) = 0; 
EUOPENSCETNDEN DS BEF 075) AND COTDCN)- GCE 0.5)); 
MXS(N)SDNO(N) = AMI(N); 
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DNO(N) = NO; 
LOOP( (N1,M)SART(N,N1,M), 
MXS(N1)SCDNOCN) + DNO(N1)) = MXS(N1) + AMT(N) 
* AMT(N1)S(NOT DN(N1)); 
DNO(N1) = ; 
1007 ) 7 777 
BGP(NI,N2,M1) — MINCBGBCNI,N2,MI); PARCONTSN2S E01 
TOBGPCNSNI,MOSCANBONIUEES0 5) FT ٦+ 77۳۲ 
MXS(N2)SCDNOCN) + DNO(N2)) = MXS(N2) -* AMT(N) 
+ AMT(N1)SCNOT DN(N1)) + AMT(N2)S(NOT DN(N2)); 
DNO(N2) = NO; 
LOOP( (N3,M2) SART(N2,N3,M2), 
BGP(N2,N3,M2) 7» MIN(BGP(N2,N3,M2), PAR(N2,N3,M2) 
* BGP(NT,N2,MI)SQMITCN2) LE 0.5) + BGT(N2 M2); 
MXS(N3)S(DNO(N) + DNO(N3)) = MAS(N3) + AMT(N) 
+ AMT(N1)S(NOT DN(N1)) + AMT(N2)S(NOT DN(N2)) 
* AMT(N3)SCNOT DN(N3)); 
DNO(N3) = NO; 
LOOP((N4,M3)SART(N3,N4,M3), 
BGP(N3,N4,M3) = MINCBGP(N3,N4,M3), PAR(N3,N4,M3) 
+ BGP(N2,N3,M2)5CAMT(N3) LE 0.5) + BGLTUNS ASIA: 
MXS(N4)S(DNO(N) + DNO(N4)) = MXS(N4) + AMT(N) 
+ AMT(N1)S(NOT DN(N1)) + AMT(N2)S(NOT DN(N2)) 
* AMT(N3)S(NOT DN(N3)) + AMT(N4)S(NOT DN(N4)); 
DNO(N4) = NO)))); 
DNO(N1)S(IND(N1) GE 0.5) = YES); 
EDP(N,N1,M)S(MN(N)S(CCAPC(M) * NUM(N,M)) GE RELT)) = MIN(EDP(N,N1,M), 
CEILCMINCMASCUN) ;MAMCN)) 7 CEAPCM) * NUMCN Mae: 
ART(N,N1,M)SCART(N,N1,M)S(CEDP(N,N1,M) - BGP(N,N1,M)) LE NRLT)) = 
ARTCN,N1,M) SCCOSUMECNZ MEJSARTONZ NAM) DIES 
AND (ABSC(AMT(N)) LE 0.5)) = NO; 
Initialize DIS, DNO and STOP if DISC equals 0. 
DIS(N,Ni)S(DISC LE 0.5) = 0; 
LOOPENMSEDISC SEE OS 
DIS(N,N1)SART(N,N1,M) = 1); 
DNO(N)S(DISC LE 0.5) 
سے‎ YESS(ABS(SMIN((N1,N2,M)S((NOT NO(N1)) * ART(N1,N2,M)), ORD(N1)) 
- ORD(N)) LE RELT); 
STOPS(DISC LE 0.5) = )نا5‎ ),۱ 01۸/١ ٤ ٣٣ 
Starting with the first nondummy node in network found, set 
DNO to 'yes' for any node which (1) is on a path with the first 
node found or (2) is on a backtrack path from other nodes on a 
path with the first node found. Stop when all arcs in the network 
have been checked. 
LOOP(IS((STOP GE 0.5) AND (DISC LE 0.5)), 
LOOP(NSDNO(N), 
LOOP(NISCCDISC(CN NI) GE ۶) )) 17705 ٣٦ 
DNO(N1) = YES; 
78ھ‎ N1)SCDIS(N, N1) GE 0.5) = 
DIS(N1,NJ)SCDIS(NI N) GE 0S o 
STOP = SUN((N, N1)SDNOCN) DISCNENI) F DIS(N1, N))); 
DNO(N)S((NOT NO(N))SCABSC(DISC - 1) LE RELT)) = YES; 
ART(N,N1,M)S(NOT (DNO(N) * DNO(N1))) = NO; 
Ensure that the first component found contains both supply and 
demand. If not, repeat the process for the next component. 
SUMASCCCSUMCNSCDNO(CN) * DN(N)), ABS(AMT(N))) LE 0.5) 
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+ + 


+ 


rij 


+ + 


OR (SUM(NS(DNO(N) * SN(N)), AMT(N)) LE 0.5))S(DISC LE 0.5)) = 
NEC(ND)SC(DNOCN)SCCSUMCNISCDNO(N1) * SN(N1)), AMT(N1)) LE 0.5) 

OR (SUUN? CUNOUCN2 so BNCN2Z)) 5. ASSCAMICN2)) LE OS) = YES; 
ARTC(N;NE,M)SCCNOT (NECCNO ^* NFCCNI))) * ARCCN,NI,N)) = YES); 


* Ensure that total supply is greater than total demand 

* in generated network. 

ABORTS((SUM(NS(DNO(N) * DN(N)), ABS(AMT(N))) 
- SUM(NSCDNO(N) ** (MN(N) * SN(N))), AHMT(N))) GE 0.5) 
"ERROR: Total demand (-) exceeds total supply (+) " 
“in generated network. ', DNO, ANT; 


* Reset the following parameters if N is not in the first connected 
* component found. 

ART(N,N1,M) = ARC(N,N1,M); 

ARC(N,N1,M)SCNOT (DNO(N) * DNO(N1))) = 

MXT(N, M)S(NOT DNO(N)) = 0; 


MXR(N,M)SCNOT DNO(N)) = 0; 
BGT(N,M)SCNOT DNO(N)) = 0; 
EDT(N,M)SCNOT DNO(N)) = 0; 


PAR(N,N1,M)S(NOT ARC(N,N1,M)) 
PRT(N,N1,M)S(NOT ARC(N,N1,M)) 
BGP(N,N1,M)S(NOT ARC(N,N1,M)) 
EDP(N,N1,M)S(NOT ARC(N,N1,M)) 
NUM(N,M)S(SUM(NISPRT(N,N1,M), 


. we we 


OOOO 


— m" Hn n 


) LE 0.5) cg m NUMCN MOS 


* Set XTE equal to 'yes' for all arcs from N to N1 vía M that 

* are part of the mine network if P is less than or equal to EDP 

* and greater than or equal to BGP. Set MNU to 'yes' if at least 

* one arc exists from N to Nl (or from N1 to N) via mode M at time 

* period P. 

ENNI PMN) = YESSC(ARC(N, NIS M)S(((ORDCP) = (EDP(N; N1, M) -+ 1)) LE RELT) 
AND CCORDCP) = BGP(N,N1,M)) GE RELT))); 

MNU(N,P)S(NOT NO(N)) = YES 
SESUHOONTSMOS dos CATE ENASNS PIM) + XTECNSNI,P + PARCN,NI,M),M))) Gê 0.5); 


* Calculate ZUP, XUP, and RHX using AMT, MXS and MXM. 
ZUP(N)SCNOT NO(N)) = ABSCAMT(N))SDN(N) * MIN(MXS(N), MXM(N))SCNOT DN(N)); 
XUP(N,N1,M)SARC(N,N1,M) = MINCZUPCN), ZUP(N1), CAP(M) * NUM(N,M)); 
RHXCN,P)S(DNO(N)SCABSCORD(P) - SMINCPISMNU(N,P1), ORD(P1))) LE RELT)) 

= ABS(AMT(N)) - 2$(ABS(AMT(N)) LE RELT); 
RHX(N,P)SCDNOCN) 

SCABS(ORD(P) - 1 = SMAX(PISMNU(N,P1), ORD(P1))) LE RELT)) = - 
RESON, F )SCDNOCN) SCRHX(N,P) GE RELI)) = MINCRHXCN,P), ZUP(N)); 


* Calculate the following parameters in preparation for forming 
* incompatible arc groups. 
MRN(N,M)SCDNO(N)SC(NUM(N,M) GE 0.5 
MEXCN,M)SCDNOCN)SCNUM(N,M) GE 0.5 
MNACN,M)SCDNOCN)SC(NUM(N,M) GE 0.5 

SCABS(PRT(N,N1,M) - MRN(N,M)) 
MNB(N,M)SCDNOCN)SCNUM(N,M) GE 0.5)) 
MXECN,M)SCDNOCN)SCNUM(N,M) GE 0.5)) SMAX(NISARC(N,N1,M), EDP(N,N1,M)); 
LSG(N,M)S$CDNO(N)SCNUM(N,M) GE 0.5)) SMAX(NISCARC(N,N1,M) 

SCCABSCPRUUN NIM) = MRN(N,M)) LE RELT) 


SMIN(NISARC(N,N1,M), PRT(N,N1,M)); 
SMAX(NISARC(N,N1,M), PRT(N,N1,M)); 
SMIN(N1$(ARC(N,N1,M) 

RELT)), PAR(N,N1,M)); 
SMIN(NISARC(N,N1,M), BGP(N,N1,M)); 
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AND (ABS(PAR(N,N1,M) - MNA(N,M)) LE RELT))), EDP(N,N1,M)); 
LSG(N,M)S(DNO(N)S(LSG(N,M) GE RELT)) = 2 + LSG(N,M) - MNB(N,M) - MRN(N,M); 
LSG(N,M)S(DNO(N)S(((MRN(N,M) - (MXE(N,M) - MNB(N,M))) GE RELT) 

AND (SUM(NISARC(N,N1,M), 1) GE 1.5))) = 1; 


* Calculate GBG and GED. This loop is explained in depth on page 
* 27 of the thésis documentation: 
LOOP((N,N1,M)S((DNO(N) * (NOT MN(N)))SC((NUM(N,M) GE 0.5) 
AND ((MRX(N,M) - 1 GE RELT) OR ((ABS(MRX(N,M) - 1) LE RELT) 
AND (SUM(N2SARC(N,N2,M), 1) GE 1.5))))), 
LOOP( PSXTE(N,N1,P,M), 
IBEG = MNB(N,M); 
LOOP(GS((LSG(N,M) - ORD(G)) GE NRLT), 
GBG(N,N1,P,M)S((GBG(N,N1,P,M) LE RELT) 
AND (ABS(PRT(N,N1,M) - MRN(N,M)) LE RELT) 
AND (ABS(PAR(N,N1,M) - MNA(N,M)) LE RELT)) 
- ORD(G)S(((ORD(P) - 1 - IBEG) GE NRLT) 
AND ((IBEG + MRN(N,M) - ORD(P) + 1) GE RELT)); 
GED(N,N1,P,M)SC((ABS(PRT(N,N1,M) - MRN(N,M)) LE RELT) 
AND (ABS(PAR(N,N1,M) - MNA(N,M)) LE RELT)) 
= MAX(ORD(G)$((ABS(IBEG + MRN(N,M) - ORD(P) + 1) GE RELT) 
AND ((ORD(P) - 1 - IBEG) GE NRLT)), GED(N,N1,P,M)); 
GBG(N,N1,P,M)S((GBG(N,N1,P,M) LE RELT) 
AND (CABS(PRT(N,N1,M) - MRN(N,M)) GE RELT) 
OR (ABS(PAR(N,N1,M) - MNA(N,M)) GE RELT))) 
- ORD(G)S(((IBEG 4 2 * MRN(N,M) - MNA(N,M) 
- ORD(P) * PRT(N,N1,M) - PAR(N,N1,M)) GE NRLT) 
AND ((ORD(P) - 1 + PAR(N,N1,M) 
- IBEG + MNA(N,M)) GE NRLT)); 
GED(N,N1,P,M)S((ABS(PRT(N,N1,M) - MRN(N,M)) GE RELT) 
OR (ABS(PAR(N,N1,M) - MNA(N,M)) GE RELT)) 
- MAX(ORD(G)S(((IBEG + 2 * MRN(N,M) - MNA(N,M) 
- ORD(P) + PRT(N,N1,M) - PAR(N,N1,M)) GE NRLT) 
AND ((ORD(P) - 1 + PAR(N,N1,M) 
- IBEG + MNA(N,M)) GE NRLT)), GED(N,N1,P,M)); 
IBEG = IBEG + 1))); 


LOOP((G,N,M)S((LSG(N,M) - ORD(G)) GE NRLT), 
LOOP((N1,P)S((GBG(N,N1,P,M) GE RELT) AND ((ABS(PRT(N,N1,M) 
- MRN(N,M)) GE RELT) OR (ABS(PAR(N,N1,M) - MNA(N,M)) GE RELT)) 
AND ((ORD(P) - PAR(N,N1,M) - MNB(N,M) - ORD(G) - MRN(N,M) 
+ MNA(N,M)) GE NRLT)), 
GBG(N,N1,P,M) = MAX(GBG(N,N1,P,M), ORD(G) + 1))); 


LOOP((G,N,M)SC(LSG(N,M) - ORD(G)) GE NRLT), 
LOOP((N1,P)S((GED(N,N1,P,M) GE RELT) AND ((ABS(PRT(N,N1,M) 
- MRN(N,M)) GE RELT) OR (ABSCPARONSNI E) = BNACNSMD) GE REGIDO 
AND ((MNB(N,M) + ORD(G) + MRN(N,M) - MNA(N,M) 
- ORD(P)- PRT(N,N1,M) * PAR(N,N1,M)) GE RELT)), 
GED(N,N1,P,M) = MIN(GED(N,N1,P,M), ORD(G) - 1))); 
ر[ہ](]()1۔۔۔۔ہ٭*٭‎ ۲)01٢ س ہہ اعد لہ بد ہبہ اعد بب اب بب اج ابع بد سس ہس ہے سے ہے ہے ہے ہے ہے ہے ہے ہے ہے سے ہے ے ے سے ہے سے ے ۔ ے ۔ سے س س س ۔ ۔ [ ]ڑاگ‎ 
POSITIVE VARIABLES 
X(N,N1,P,M) flow of mines from N to Nl via mode M at time period P 
Z(N,P) mines inventory at N at time period P; 
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* fix upper bounds on positive variables; set some positive inventory 
* variables to be constant 

SIUPONSNISP.HDSXTECONONISPOI) SS RUPONSNI SO): 

. UP(N,P) = ZUP(N); 


N 


FREE VARIABLE 
۷ objective value; 


BINARY VARIABLE 
۰۲۱۰۰۱۷۰٠۷۶۱٦۱۰۰ ۱1۱۱071176۲0 6.017: lnot )- used from N to NI-via Mat P; 


EQUATIONS 
OBJV objective value 
MNFB(N,P) flow balance equations for mine network 
LINK(N,N1,P,M) link equations between mine and mode networks 
MDSK(N,M,G) mode schedule equations for incompatible arcs 
MDPG(N,M) constraint to limit total time that M is from N 
MDRG(N,M) constraint to limit total runs that M is from N; 
OB IV 2% 


V =G= SUM((N,N1,P,M)SXTEC(N,N1,P,M), CPRN(N) * PRN(N1) * PRM(M) 
O OORT SI TTR NJ 2 T PON NISSEN); 


P)S(MNU(N,P) * DNO(N)).. 

O =G= (1SDN(N) - 1$(NOT DN(N))) 

* (Z(N,P)SCCREX(N,P) LE 0.5) AND (RHX(N,P) GE -1.5)) 

- Z(N,P + 1)5CREX(N,P + 1) GE -0.5) * RHX(N,P)S(RHX(N,P) GE 0. 5)) 
SUMC(N1,M), X(N1,N,P,M)SXTE(N1,N,P,M) 
- X(N,N1,P + PAR(N,N1,M),M)SXTE(N,N1,P + PAR(N,N1,M),M)); 


DEEST ONT PSMYSATEEN SNIT PM Jay 
DESDE AUPCNOBAEMQ C YUNSNIZPIMJ -OAUONSNISP;MD; 


MDSK(N,M,G)SCDNOCN)SCCLSG(N,M) - ORD(G)) GE NRLT)).. 
] 2G- SUM((N1,P)SCCCORD(G) - GBG(N,N1,P,M)) GE NRLT) 
ABDOSCCGEDCNSNT;PSM) - ORD(CG)O GE NRET)), YCON,NIS,P M5); 


HUPBSCNSMOSCCDNOUNO Se CNOT CNOCN) * MOCHOO))SCCTEND. -*- HMXICNSM)O GE RELT)).. 
SEMCONISPOSNXIECNAINISPSMDSSPRTCN,NTSM) * YON NIGP NJ) ZI MATCN,M)S 


DOSOONSMISCCBNOCNO OS (NOT CNOCN) * MOCM))))SCMARCN,M) GE RELT)).. 
BEUCUNISE)SXTECNSNISPSMOS, YCN,NJ4PSM)) —L-9 MARCNSMOS 


MODEL SUMIT /OBJV,MNFB,MDSK,LINK,MDPG,MDRG/; 


SOLVE SUMIT USING MIP MINIMIZING V; 


D O N cl lec ا یک‎ e ES DERE 
DISPLAY " The following arcs were considered by SUMIT.", ARC; 

DISPLAY " *** All times in given in hours, incremented by TPER. ***"'; 
PARAMETER 


TIMECOMP(N, X) Time in which layin, of mine field N was completed; 
TIMECOMP(N, ' TIME. DONE' کت‎ DNO(N)) 
= SMAX((N1, DOMOSCOYSDONI NSTPIM) GESRBBEIO. TBER C CORDCP)-*. 1)95; 


95 


OPTION TIME COUR. 
DISPLAY ٦ 


PARAMETER 
MOMAGREP(N,*) MOMAG Schedule Report (Start Time is 0); 
MOMAGREP(N, START. TIME’ )S((MN(N) * DNO(N)) 
SCSMAXCCNI,P,M)SXTECN,NI, P.M). UNNI EDS GE TEND 
= TPER * SMINCCNI, P MSC. LCN, NISP M) GE RELT), ORDCD) = 2): 
MOMAGREP(N,' STOP. TIME' ) CCHNCN) * DNO(N)) 
S( SMAX((NI, P,M)SXTECN,N1, ۶٘ٗ IDD CE ۱ ٦١ 
= TPER * SMAX((N1, P :M)SCX. LN, NIE) GE RELT), ORD(P) - 1); 
MOHAGREP(N, ' NUM. MINES' )SCMNCN) * DNOCN)) 
- SUMC(CNISPS;MOSCX. L(NSNISPSSDSGESRESDTOSEASDONSONTSPS MS 
OPTION MOMAGREP: 1: 1: 1; 
DISPLAY MOMAGREP; 


PARAMETER 
SCHEDULE(N,M,G,N1,*) Schedule for Mine Transshipment by Origin Base; 
LOOP(NS((DNO(N) * (NOT MN(N))) 
S(SUM((N1,P,M), X. L(N,N1,P,M)) GE RELT)), 
LOOP(MS(NUM(N,M) GE 0.5), 
IBEG = 1; 
LOOP((P,N1)S(X. L(N,N1,P,M) GE RELT), 
LOOP(GS(ABS(ORD(G) - IBEG) LE RELT), 
SCHEDULE(N,M,G,N1, 'NO. MODES' ) 
= CEIL(X. L(N,N1,P,M) / CAP(M)); 
SCHEDULE(N,M,G,N1, , NO. MINES") - X. L(N,N1,P,M); 
SCHEDULE(N,M,G,N1, ' LOAD. TIME') 
= TPER * (ORD(P) - 1 - PAR(N,N1,M)); 
SCHEDULE(N,M,G,N1, ETD. ORIG') 
= (TPER * ( ORD( P) - 1 - PAR(N,N1,M))) + TEV(M,'TLD' ); 
SCHEDULE(N,M,G,N1,'ETA. DEST' ) 
= (TPER * (ORD(P) - 1)) - TEV(M,'TUL’); 
SCHEDULE(N,M,G,N1, AVAIL. TIME’ ) 
= TPER * (ORD(P) - 1 + PRT(N,N1,M) - PAR(N,N1,M))); 
IBEG = IBEG + 1)); 
OPTION SCHEDULE: 1:4: 1; 
DISPLAY SCHEDULE; 
SCHEDULE(N,M,G,N1, ' NO. MODES') = 
SCHEDULE(N,M,G,N1, ' NO. MINES' NS 
SCHEDULE(N,M,G,N1,' LOAD. TIME ' ) = 
SCHEDULE(N,M,G, N1, 'ETD. ORIG') = A 
SCHEDULE(N,M,G,N1, ES DEST JD O 
SCHEDULE(N,M,G,N1, AVAIL. TIME’) = 


PARAMETER 
MINEINV(N,*) The number of mines in left in inventory; 
MINEINV(N, ‘MINES. LEFT' )$(DNO(N) * (MN(N) + SN(N))) 
= SUM((N1,P,M)S((X. L(N,N1,P,M) GE 0.5) OR (X. L(N1,N,P,M) 4 ۵0٣ 
K. LONIGN,P,M) = KX LNSNI,P. NIJ FT ANI 
OPTION MINEINV: 0; 
DISPLAY MINEINV; 


SET 


NEXTRUN(N,*) Nodes that should be included in the next run of SUMIT; 
NEXTRUN(N, ' NEXT. RUN' )SC(NOT (NO(N) + DNO(N) + NFC(N))) 


56 


$((SMAX((N1,M) SART(N,N1,M), ABS(NUM(N,M))) GE 0.5) 
OR (AMT(N) LE -0.5))) = YES; 
LOOP(IS((ABS(ORD(I) - 1) LE 0.5) 
AND ( SUH(NSCNEXTRUN(N , " NEXT. RUN' ) * DN(N)), AMT(N)) LE -0.5) 
AND (SUM(NSNEXTRUN(N, NEXT. RUN'), 1) GE 0.5)), 
DISPLAY " Change NUM(N,M) and DN(N) to O for all nodes not in", 
'" NEXTRUN(N) and rerun SUMIT to finish the solution.", NEXTRUN; 
LOOR( TISCCABSCORDC TI) 909) LE 0:5) 
AND (SUM(NSNEXTRUN(N, NEXT. RUN), AMT(N)) LE -0.5)), 
AMT(N) 2» AMT(N)SNEXTRUN(N, ' NEXT. RUN' ); 
OPTION AMT: 0; 
DISPLAY "WARNING: There is not enough supply to meet demand ", 
"for the next run of SUMIT. ", AMT)); 
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APPENDIX B. SAMPLE OUTPUT, FIRST RUN 
COMPILATION TIME = 0.950 SECONDS 


MODEL STATISTICS 


BLOCKS OF EQUATIONS 6 SINGLE EQUATIONS 173 
BLOCKS OF VARIABLES > SINGLE VARIABLES 213 
NON ZERO ELEMENTS 625 DISCRETE VARIABLES 78 


GENERATION TIME 1.190 SECONDS 


EXECUTION TIME 


Ii 
SI 


.500 SECONDS 


SOLVE SUNNET 
MODEL SUMIT OBJECTIVE y 
TNBE MIP DIRECTION MININIZE 
SOLVER ZOOM FROM LINE 821 
Y SOLVER STATUS 1 NORMAL COMPLETION 
Wii MODEL STATUS 1 OPTIMAL 
sare OBJECTIVE RVA BUE 45612.0000 
RESOURCE USAGE. AUT 158. 146 10000. 000 
ITERATION COUNT, LIMIT 763589 10000000 
ZO Ot یس‎ ٣ --- Version 2.1 Jun 1955 


Courtesy of Dr Roy E. Marsten, 
Department of Management Information Systems, 
University of Arizona, 
Tucson Arizona 85/21; U.S.A 


No options file found - using defaults. 


Work space needed (estimate) -- 67966 words. 
Work space available -- 67966 words. 
Maximum obtainable -- 294142 words. 


The LU factors occupied 804 slots (estimate 3068). 
The branch and bound tree contained 4550 nodes (max. 10000 nodes). 


Iterations: Initial LP 193, Time: 0.26 


5S 


ARCS IN NETWORK FROM NODE N TO NODE 
Ni VIA MODE M 


3214, Beo] 

73061, 149.45 

120 0.14 
NONOPT 


P3 B52 TACAIR SSN688DM DBRM 


YES 


YES 


YES 


YES 


YES 
YES 


YES 
YES 
YES 


YES 
YES 


Heuristic 
Branch and bound 
Final LP 


mera REPORT SUMMARY : 0 
O INFEASIBLE 
O UNBOUNDED 


THE FOLLOWING ARCS WERE CONSIDERED BY SUMIT. 


---- 815 


Eme S15. SEI ARC 


TRUCK C130 


YES 
NES 


. S8GQ 

. ANDERSNAFB YES 
. S9RP 

. MCASIWKNI 

. MINEFIELD3 

. $10JA 

. MINEFIELD2 

. MINEFIELD3 


NASADAK 
MCASIWKNI 
MINEFIELD] 
MINEFTELDZ 
MINEFIELD3 


. 71 
.. 2 
+م,‎ 27 
SMINEFTELDI 


rxe ALL TIMES IN GIVEN-IN HOURS, INCREMENTED BY- TPER. rr 


TIME IN WHICH LAYING OF MINE FIELD 
N WAS COMPLETED 


MOMAG SCHEDULE REPORT 
(START TIME IS 10) 


M8GQ 
S86Q 
M9RP 
S9RP 
S9RP 
M10JA 
S10JA 
S10JA 
ANDERSNAFB. 
ANDERSNAFB. 
ANDERSNAFB. 
ANDERSNAFB. 
ANDERSNAFB. 
NASADAR 
MCASIWKNI 
CVNA 
SSN688A 


---- 817 


---- 824 PARAMETER TIMECOMP 


TIME. DONE 


36:0 
36.0 
56. 0 


MINEFIELD} 
HISEFIELD2 
MINEFIELD3 


---- 837 PARAMETER MOMAGREP 


STOR TIKE NUM. MINES 


SCHEDULE FOR MINE TRANSSHIPMENT 


BY ORIGIN BASE 


LOAD. TIME ETD.ORIG ETA.DEST AVAIL. TIME 


6.0 
120 
18.0 
24.0 


TIME 


18.0 
30-0 
42.0 


4.0 
10.0 
16.0 
2 20 


1.0 
6.0 720 
1250 1380 
1820 1950 


SCHEDULE FOR MINE TRANSSHIPMENT 


BY ORIGIN BASE 


LOAD. TIME ETD.ORIG ETA.DEST AVAIL. 


DT 
257 
3537 


6.0 7 
18.0 19.0 
30.0 31 0 


SCHEDULE FOR MINE TRANSSHIPMENT 
BY ORIGIN BASE 


LOAD. TIME ETD. ORIG ETA.DEST AVAIL. TIME 


1.0 2:9 18.0 
18.0 19. 0 2357 3050 
3080 31550 321 42.0 


SCHEDULE FOR MINE TRANSSHIPMENT 
BY ORIGIN BASE 


LOAD. TIME ETD. ORIG ETA. DEST 00 TS 


4.0 11.5 240 


24.0 28.0 DOD 48.0 


SCHEDULE FOR MINE TRANSSHIPMENT 
BY ORIGIN BASE 


LOAD. TIME ETD. ORIG ETA.DEST AVAIL. TIME 


120 2 720 6.0 
18.0 7 7ا 13.0 12.0 
24.0 7 28% 19519 18.0 


60 


54.0 
1250) 
54.0 


860 PARAMETER SCHEDULE 


NO. MINES 
20 
20.0 


20.0 
20.0 


SCHEDULE 


NO. MINES 
48.0 


48.0 
3280 


SCHEDULE 


NO. MINES 
48.0 


48.0 
48.0 


SCHEDULE 


NO. MINES 
۵ھ‎ 
145.0 


SCHEDULE 


NO. MINES 


SEIO 
51.0 
51.0 


M8GQ 650 
M9RP 1239 
M10JA 6-0 
S8GQ. TRUCK 

NO. MODES 
1. ANDERSNAFB 1.0 
2. ANDERSNAFB 0 
3. ANDERSNAFB 1:0 
4. ANDERSNAFB DU 
ea 850 PARAMETER 
SIRP. P3 

NO. MODES 
1. MINEFIELD3 6.0 
2. MINEFTELD3 60 
3. MINEFIELD3 4.0 
reme 860 PARAMETER 
S10JA PS 

NO. MODES 
3-4۲ 609 
2. MINEFIELD2 6.0 
3. MINEFIELD2 6.0 
pcm 860 PARAMETER 
ANDERSNAFB. B52 

NO. MODES 
1. MINEFIELDS 20 
Ze TENS ۵ 70 
ےج میم‎ 860 PARAMETER 
CVNA. TACAIR 

NO. MODES 
8 77 150 
2. MINETIELDZ 150 
 ؛‎ 5 +4 1.0 


se 860 PARAMETER SCHEDULE SCHEDULE FOR MINE TRANSSHIPMENT 
BY ORIGIN BASE 


SSN688A. SSN688DM 
NO. MODES NO.MINES LOAD. TIME  ETD.ORIG ETA.DEST AVAIL. TIME 


1۰۰۲٣٢ 10 50.0 12:9 12.0 225 30.0 


E 874 PARAMETER MINEINV THE NUMBER OF MINES IN LEFT 
IN INVENTORY 


۹۷۷٣ 


M8GQ 146 
M9RP 1 
M10JA 216 
S10JA s 
CVNA 8 


ee. 883 CHANGE NUM(N,M) AND DN(N) TO O FOR ALL NODES NOT IN 
NEXTRUN(N) AND RERUN SUMIT TO FINISH THE SOLUTION. 


eee D53 SET NEXTRUN NODES THAT SHOULD BE INCLUDED 
IN THE NEXT RUN OF SUMIT 

NEXT. RUN 

M11SC YES 

51150 YES 

NASJAX YES 

CVNB YES 

SSN688B YES 

MINEFIELD4 YES 


xs FILE SUMMARY FOR USER 8693P 


INPUT TESTI GAMS Å 
OUTPUT TESTA LISTING A 


EXECUTION TIME = 1. 370 SECONDS 
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APPENDIX C. SAMPLE PROGRAM, SECOND RUN 


STITLE * * * Scheduler for Mine Transshipment (SUMIT) * * * 

*- ---DOCUMENTATION===-----=------- oca canenannanan aaa nes 
Software to Optimize Schedules for Mine Transshipment 

Developed for Commander, Mine Warfare Command 

Written by Tammy L. Glaser, LT, USN 

Complete documentation contained in thesis 

06 September 1991 

*----GAMS AND DOLLARS CONTROL OPTIONS----------------------2--------2-2-2222- 
* Do not change the following options. 

SOFFUPPER OFFSYMXREF OFFSYMLIST 

OPTIONS SOLPRINT = OFF, LIMROW = 0, LIMCOL = 0; 

* Increase the following options only when recommended in the 

* solution report of the output. 

OPTIONS RESLIM = 10000, ITERLIM = 10000000, WORK = 10000; 

* Set OPTCR as recommended on page 72 of the thesis documentation. 
OPTIONS OPTCR = 0.00001; 


voco ook 


*----DEFINITIONS AND DATAÃ---==================smenseene eee... O 
* See Chapter 2 of thesis documentation to explain data in detail. 
SCALARS 
* The following scalars must be set to values greater than 0. 
TEND end time of problem (h) [2299 
TPER length of time period (h) 7 
RELT zero tolerance for comparisons of real data /0.001/ 
* The following scalars are optional (enter 0 for default value). 
MNLA minimum distance for air modes over land routes (nm) /50/ 


MXSS maximum distance for transfers between supply nodes (nm) /100/ 
DISC 1 to run entire problem and 0 to split into subproblems  /0/; 


SETS 
The first two sets must have the element O in them. Note that 

N must contain all elements in CN(N), found in the optional section. 

M must contain all elements in DM(M) (and TM(M), BM(M), CM(M), and 

LM(M), which are found in the optional section). See page 7 of 

thesis documentation to explain the MOMAG-supply node split. 

N nodes /0, M8GQ, S5CQ, M9RP, S9RP, MIOJA; ۳ء یھ‎ 
M11SC, S11SC, ANDERSNAFB, NASADAK, 
MCASIWKNI, NASJAX, CVNA, CVNB, SSN688A, 
SSN688B, MINEFIELD1, MINEFIELD2, 
MINEFIELD3, MINEFIELD4/ 


$$ + + + 


M modes /0, TRUCK, C130, C141, 23, 8 ۰ 
SSN688DM, DBRM/ 
DM(M) delivery modes /P3, B52, TACAIR, SSN688DM/ 


* P must start with O and end with the value of the number of time 
periods plus one. G and I must start with 1. G must end with 
* the number of time periods, while I must end with the square of 


+ 


* the number of nodes (nonzero elements in N). 
P time periods ug 9 235 
G group numbers nct 22 
I iteration numbers /1 * 400/ 


* The following sets are optional. Enter 0 between slashes for sets 
* not needed for the problem. If CN(N) is empty (i.e., initialized to 


* 0), nodes are assumed to be land nodes. If there are no transshipment 
* modes or disassembled mines, enter O between slashes for TM(M) or 

* BM(M). A11 modes are assumed to be air modes if CM(M) and LM(M) 

* are initialized to O. 


CN(N) sea nodes /CVNA, CVNB, SSN688A, SSN688B/ 
TM(M) transportation modes T TRUCK + +80 77 
BM(M) dummy build rate modes / DBRM/ 
CM(M) dummy sea modes / SSN688DM/ 
LM(M) land modes [TRUCK / 
* Do not change the following evolution set. 
E mode evolution FI EDS TULL, TODA: 


* Do not change the following set aliases: 
ALIAS (N,N1,N2,N3,N4); 

ADDAS (M,NM1,M2,M3); 

ALIAS (P,P1); 

ALIAS (I TT); 


PARAMETERS 
* The following parameter is required for node N: 
AMT(N) supply (+) or demand (-) at node N (mines) 
/M8GQ 200 
S8GQ 25 
MSRP 159 
M10JA 270 
S10JA 95 
M11SC 100 
5115C 40 
ANDERSNAFB 89 
CVNA 161 
CVNB 100 
SSN688A 0 
SSN688B 40 
MINEFIELDI 0 
MINEFIELD2 0 
MINEFIELD3 0 
MINEFIELD4 -201/ 


* The following parameters are optional (leave blank line with no 
* slashes for default values): 


PRN(N) priority of node N 
ZINSEIELDI 2 
MINEFIELD? 1 
MINBSEIBEDS 3 
۶ 1/ 
MXD(N) maximum distance between mine fields at node N (nm) 
/MINEFIELD2 50/ 
BGD(N) transit distance for sea node N at beginning (nm) 
/CVNA 33 
CVNB 75/ 
EDD(N) transit distance for sea node N at end (nm) 
/CVNB 75/ 
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transit speed for sea node N (nm per h) 


25 
25/ 


/CVNA 
CVNB 


SPXCN) 


The following parameters are required for mode M: 
capacity of mode unit in mode group M (mines per mode) 


/ TRUCK 20 
C130 20 
CIS 20 
P3 8 
B52 22 
TACAIR 51 
SSN688DM 50/ 

maximum distance that mode M may travel (nm) 

/ TRUCK 150 
C130 3000 
C141 3000 
P3 1500 
B52 3000 
TACAIR 500 
SSN688DM 200/ 

transit speed of mode M (nm per h) 

/ TRUCK 30 
C130 380 
C141 410 
P3 350 
B52 480 
TACAIR 450 
SSN688DM 15 
DBRM 9/ 


CAP(M) 


MXL(M) 


SPD(M) 


The following parameter is optional (leave blank line for default): 


C C) I5‏ مم مم مم UD‏ دم 


/ 


priority of mode M 
/ TRUCK 

C130 

C141 

P3 

B52 

TACAIR 

SSN688DM 

DBRM 


PRM(M) 


The following parameter is required for mode M at node N: 
number of mode M at node N (modes) 


FUO0rO0Oo0oOo0oOo0oooOo 


/M8GQ. DBRM 
S8GQ. TRUCK 
M9RP. DBRM 
S9RP. C130 
S9RP. P3 
M10JA. DBRM 
S10JA. P3 
M11SC. DBRM 
S11SC. TRUCK 
ا‎ 


NUM(N,M) 


O O ©‏ نین یی تن بتپرں پیم 


/ 


ANDERSNAFB. B52 
ANDERSNAFB. C130 
۸۸3 
MCASIWKNI. TACAIR 
NASJAX. P3 

CVNA. TACAIR 
CVNB. TACAIR 
SSN688A. SSN688DM 
SSN688B. SSN688DM 


* The following parameters are optional (leave blank line for default): 
earliest beginning time for mode M at node N (h) 
/ SSN688A. SSN688DM 


10/ 


latest end time for mode M at node N (h) 
/ CVNA. TACAIR 


24/ 


maximum total time mode M can be absent from node N (h) 


maximum number of runs for mode M from node N 
/ ANDERSNAFB. B52 
ANDERSNAFB. C130 
SIRE, C130 


2 
2/ 


1 means eliminate runs from nodes N to N1 via mode M 
T99RP.MINEFIELD2. F3 1/ 


S10JA 
3(٣ 
1157 

738 
738 
20:5 
0 
7028 
7028 
1122 
2561 
1083 
1032 
849 
2293 
2144 
/388 
2170 
1015 
1203 
7414 


M10JA 
1157 
11597 

738 
78 
0 
EU 
7028 
7028 
1172 
2881 
1083 
1032 
849 
1295 
2144 
7388 
2170 
1015 
1203 
7414 


BGT(N,M) 


EDT(N,M) 


MAIN.) 


MAR(N,M) 


)01 + 8ھ 


* The following two parameters are required. 
* READ PAGE 22 OF THESIS DOCUMENTATION, WHICH EXPLAINS WHEN NEGATIVE 
ZODISTANCES MUST BE USED. 
TABLE DIS(N,N1) distance between node N and node N1 (nm) 


S9RP 
1345 
1345 
05 
0 
738 
738 
1754 
7754 
1364 
3614 
1521 
7768 
1580 
8020 
2881 
5125 
2909 
1725 
897 
5.157 


M9RP 
1345 
1345 
0 
UNS 
TIB 
738 
7754 
7754 
1364 
3614 
1821 
7768 
1580 
8020 
2881 
8125 
2909 
1725 
897 
8152 


S66Q 
2055 
0 
1345 
1345 
197 
00 
7501 
7301 
>0 
2059 
1651 
7485 
1566 
7789 
2405 
7612 
2459 
1819 
2207 
27 


M8GQ 
0 


2055 
1345 
1345 
Lis, 
1157 
1501 
DOT 

=20 
2939 
1631 
7485 
1566 
1789 
2405 
7812 
2459 
1819 
2207 
42937 


M11SC 5115 ANDERSNAFB NASADAK MCASIWKNI 


1031 


2359 


=30 


7501 


7 


M8GQ 

S8GQ 

M9RP 

S9RP 

M10JA 
S10JA 
M11SC 
S11SC 
ANDERSNAFB 
NASADAK 
MCASIWKNI 
NASJAX 
CVNA 

CVNB 
SSN688A 

6 632 
MINEFIELD1 
MINEFIELD2 
MINEFIELD3 
MINEFIELD4 


4 
M8GQ 


=30 
1364 
1364 
1172 
2 
1520 
7520 
0 
2939 
1653 
7500 
1583 
7801 
2421 
7831 
2455 
1800 
2188 
7812 


2939 1631 
3614 1821 
3614 1821 
2881 1083 
2881 1083 
4222 5942 
4222 5942 
2939 1653 
0 1824 
1824 0 
4209 5950 
2088 267 
4511 5942 
777 1067 
4546 6305 
744 2059 
2032 SD 
3850 2053 
4551 6305 
SSN688A  SSN688B 
2405 2ھ‎ 
2405 7812 
2881 8152 
2881 8152 
2144 7414 
2144 7414 
4899 420 
4899 420 
2421 7831 
777 4546 
1067 6305 
4900 366 
1324 6546 
5164 327 
0 5249 
5249 0 
110 5220 
1244 6413 
3073 8023 
5270 174 
MINEFIELD3 MINEFIELD4 
2207 77 
2207 7797 
897 5152 
897 8152 
1205 7414 
1203 7414 
7607 548 
7607 548 
2188 7812 
3850 4551 
2033 6332 
7660 463 
1786 6578 


CVNB 
7789 
7789 
8020 
8020 
7293 
7.293 
294 
294 
7801 
7801 
6216 
351 
6445 
0 
5164 
327 
5146 
6294 
1817 
500 


MINEFIELD2 
1819 
1819 
1725 
18725 
1015 
TOTS 
6032 
0032 
1800 
2032 

352 
6051 
234 


66 


2501 7501 
1754 7754 
1754 1754 
7028 1026 
7028 7028 
0 20.5 
2085 0 
7520 7520 
4222 4222 
42 5942 
-124 -124 
6175 6175 
294 294 
4899 4899 
420 420 
4866 4866 
60327 6032 
7607 27607 
548 548 
NASJAX ۵۸۵ 
7485 1566 
7485 1566 
7768 1580 
7768 1580 
2092 849 
7032 849 
-124 6175 
-124 6175 
7500 1583 
4209 2088 
5950 267 
0 6188 
6188 0 
366 6445 
4900 1324 
366 6546 
4869 1344 
6051 234 
7660 1786 
463 6578 
MINEFIELD1 
2459 
2459 
2909 
2909 
2170 
2170 
4866 
4866 
2455 
744 
1090 
4869 
1344 


S86Q 

M9RP 

SIRE 

MIOJA 
S10JA 
M11SC 
51150 
ANDERSNAFB 
NASADAK 
100 +1 
NASJAX 
CVNA 

CVNB 
SSN688A 
SSN688B 
MINEFIELD1 
MINEFIELD2 
MINEFIELD3 
MINEFIELD4 


+ 
1860 
80 
M9RP 
S9RP 
M10JA 
S10JA 
M11SC 
S115€ 
ANDERSNAFB 
NASADAK 
MCASIWKNI 
NASJAX 
CVNA 
CVNB 
SSN688A 
SSN688B 
MINEFIELD1 
MINEFTELD2 
MINEFIELDS 
MINEFIELD4 


+ 
M8GQ 

5860 

MORE 

SIRE 

M10JA 
5104 
M115C 
S115C 
ANDERSNAFB 
NASADAK 
MCASIWKNI 
NASJAX 
CVNA 


CVNB 5146 6294 7817 500 
SSN688A 110 1244 3073 5270 
SSN688B 5220 6413 8023 174 
MINEFIELDI 0 1265 3081 5243 
HINEFTELD2 1263 0 1817 6459 
MINEFTEED3 3081 1817 0 8109 
MINEFIELD4 5243 6459 8109 82 
TABLE TEV(M,E) time needed to complete evolution E for mode M (h) 

TLD TUL TGD 
TRUCK 1.00 2. 00 1-00 
C130 4.00 2:00 8. 00 
C141 4.00 2.00 8. 00 
p3 1.00 0. 25 4. 00 
B52 4.00 00 8. 00 
TACAIR 1:00 0525 2. 00 
SSN688DM 0.00 1:50 0. 00; 
yere esee ze seges eve eese eese e deve dese eee ec desee de vede dee de dede oe ee ve dese dee deve ee ye dede dec Je de de eee es ede dede de de de ye eye des 
= y 
* Input is complete. Do not change any statements after this line. * 
* * 


dese dede de dese dede de rede de de dee dese de‏ جاو باو جلو لو eye eyes cv deve yes ey eyes eve eee desc Jes eye deve de Je Je Jes Jede Je Jew e Je Jedes ye Je Je Je‏ پا او از مج 
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APPENDIX D. SAMPLE OUTPUT, SECOND RUN 
COMPILATION TIME = 0.950 SECONDS 


MODEL STATISTICS 


BLOCKS OF EQUATIONS 6 SINGLE EQUATIONS 123 
BLOCKS OF VARIABLES 4 SINGLE VARIABLES 154 
NON ZERO ELEMENTS 446 DISCRETE VARIABLES 51 


GENERATION TIME 1. 170 SECONDS 


EXECUTION TIME 7.440 SECONDS 


SOLVE SUMMARY 
MODEL SUMIT OBJECTIVE SIY 
TYPE MTE DIRECTIONMININIMIZE 
SOLVER ZOOM FROM LINE 821 
٭٭٭٭‎ SOLVER STATUS 1 NORMAL COMPLETION 
sed MODEL STATUS 1 OPTIMAL 
veio OBJECTIVE VALUE 158139. 0000 
RESOURCE USAGE, LIMIT 1055820 10000. 000 
ITERATION COUNT TEINIT 66392 10000000 
"ھ2‎ ۶۹ AX MOL 2 Version 2.1 Jun 1988 


Courtesy of Dr Roy E. Marsten, 
Department of Management Information Systems, 
University of Arizona, 
Tucson Arizona 85021. 57 


No options file found - using defaults. 


Work space needed (estimate) -- 62593 words. 
Work space available == 62593 words: 
Maximum obtainable -- 246526 words. 


The LU factors occupied 552 slots (estimate 2101). 
The branch and bound tree contained 4106 nodes (max. 10000 nodes). 


Iterations: Initial LP 118, Time: OFT2 
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Heuristic 2.53.5 0-562 


Branch and bound 65947, 104.94 
Final LP 74, 0. 07 
RR REPORT SUMMARY : 0 NONOPT 


0 INFEASIBLE 
O UNBOUNDED 


815 THE FOLLOWING ARCS WERE CONSIDERED BY SUMIT. 


815 SET ARC ARCS IN NETWORK FROM NODE N TO 


NODE N1 VIA MODE M 


TRUCK C130 P3 TACAIR SSN688DM DBRM 


"c SIISC YES 
S11SC .NASJAX YES YES 
NASJAX . MINEFIELD4 YES 
CVNB .MINEFIELD4 YES 
SSN688B. MINEFIELD YES 
---- 817 *** ALL TIMES IN GIVEN IN HOURS, INCREMENTED BY TPER. *%% 
~--- 824 PARAMETER TIMECOMP TIME IN WHICH LAYING OF MINE FIELD N 
WAS COMPLETED 
TIME. DONE 
MINEFIELD4 19.2 
---- 837 PARAMETER MOMAGREP MOMAG SCHEDULE REPORT 
(START TIME IS 0) 
STOP. TIME NUM. MINES 
M11SC 3. 6 21.0 
---- 860 PARAMETER SCHEDULE SCHEDULE FOR MINE TRANSSHIPMENT 


BY ORIGIN BASE 


S11SC. TRUCR 


NO. MODES NO. MINES LOAD. TIME  ETD.ORIG ETA.DEST AVAIL. TIME 


20 40.0 150 2و5‎ 1352 
200 21:50 36 4.6 8.8 16. 6 
860 PARAMETER SCHEDULE SCHEDULE FOR MINE TRANSSHIPMENT 
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1. NASJAX 
2. NASJAX 


BY ORIGIN BASE 


) P3 

NO. MODES NO. MINES LOAD. TIME ETD.ORIG ETA. DEST AVAIL. TIME 
1. MINEFIELD4 5.0 40.0 2 8. 2 0,0: 1556 
2. MINEFIELD4 3.0 210 15.6 16. 6 1859 24.0 
s 860 PARAMETER SCHEDULE SCHEDULE FOR MINE TRANSSHIPMENT 


BY ORIGIN BASE 
CVNB. TACAIR 
NO. MODES NO. MINES LOAD. TIME ETD. ORIG ETA. DEST AVAIL. TIME 


1. MINEFIELD4 1.0 5170 3.6 4.6 7 9506 
2. MINEFIELDS 1.0 49.0 96 10.6 Ie 3 
داد ھن‎ 860 PARAMETER SCHEDULE SCHEDULE FOR MINE TRANSSHIPMENT 


BY ORIGIN BASE 
SSN688B. SSN688DM 
NO. MODES | NO. MINES ETA.DEST AVAIL. TIME 


1. MINEFIELD4 1.0 40. 0 Ip 7 25d 
= 874 PARAMETER MINEINV THE NUMBER OF MINES IN LEFT 
IN INVENTORY 
MINES. LEFT 
11150 79 


we FILE SUMMARY FOR USER ۶۳ 


INPUT TEST? GAMS À 
OUTPUT TEST? LISTING A 


EXECUTION TIME = 1. 140 SECONDS 
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APPENDIX E. GAMS IMPLEMENTATION OF SUMIT 


A. SUMIT OUTPUT 

The output of SUMIT is composed of standard GAMS reports and of data dis- 
played by the program. The standard GAMS reports include an echo print of the pro- 
gram, compilation time, model statistics (1.c., number of constraints and continuous and 
discrete variables), and solve summary. The echo print of the program was deleted from 
Appendix B since the program is given in Appendix A. The data output by SUMIT was 
also edited to make Appendix B fit within the page margins. The data displayed by the 
program include the arcs generated for the mine network (which has the same arcs as the 
mode network), the time in which the mine field deployment was complete, the 
MOMAG schedule report, the schedules for mine transshipment, the mine inventory at 
the end of the problem, and the nodes that should be included in the next run of 
SUMIT. SUMIT prints out the mine network to allow the user to verify that all arcs 
in the mine and mode networks are appropriate. The MOMAG schedule report con- 
tains the times in which each MOMAG must start and stop assembling mines and the 
number of mines that each MOMAG must build. SUMIT then prints out separate mine 
transshipment schedules for each origin node that reports the following information: 
Origin node, mode group, run number for the node-mode combination, destination node, 
total number of mines carried by the mode group, number of mode units needed, the 
time to start loading the mode group, estimated time of departure from the origin node, 
estimated time of arrival at the destination node, and the time in which the mode is 
available for loading for its next flight. It also warns the user if total supply is less than 


total demand for nodes recommended in the next run. 


B. END TIME AND TIME PERIOD LENGTH 

To find appropriate values for the end time and time period length for a new prob- 
lem, the user should take the following steps: the user should initially execute GAMS 
using the relaxed version of ZOOM, which is faster since it allows binary variables as- 
sume anv Values between zero and one. To indicate the relaxed version of ZOOM, the 
user must change the problem identifier “MIP” to “RMIP” in the “SOLVE” statement, 
which is found on line 821 of the program (see Appendix A, page 55). The goal of 
making these quick runs is to identify the latest end time with an appropriate number 


of periods in which the solver reports an infeasible solution (“INFEASIBLE”) in the 
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section of the output labeled "SOLUTION REPORT." Then, the user should run 
SUMIT using an end time equal to the latest infeasible end time plus one extra time 
period with binary requirements enforced (replacing "RMIP" with “MIP” in the 
"SOLVE" statement). In most cases, adding an extra time period vill make the solution 
feasible. If the problem is not feasible when the binary variables are enforced, the 
“SOLUTION REPORT” portion of the output will contain the phrase “INTEGER 
INFEASIBLE”. If the solution is integer infeasible, the user should first verify that 
there is enough supply to meet the demand for the mine network output by the SUMIT. 
If the supply is sufficient, then the user should rerun the “MIP” version of ZOOM after 
increasing the end time by one extra time period until a feasible solution is found. Since 
the purpose of these final runs is to find the first integer feasible solution, the user may 
set the OPTCR at 1.0 to find the solution as quickly as possible. 


C. OBTAINING NEARLY OPTIMAL SOLUTIONS 

A general rule of thumb is recommended to obtain nearly optimal solutions for 
SUMIT. When the user makes the first run of the “MIP” version of ZOOM, the user 
should set OPTCR equal to 0.10 in line 16 of the program, rerun SUMIT, and verify that 
the reported solution is satisfactory. If the number of iterations 1s less than 10,000, 
which is given in the “SOLUTION REPORT” portion of the output indicated by “IT- 
ERATION COUNT”, rerunning SUMIT with a very small OPTCR, such as 0.01, will 
probably vield an optimal solution within a reasonable number of iterations, but will not 
usually prove it. If the number of iterations is greater than 10,000, the user may rerun 
SUMIT at an OPTCR between 0.01 and 0.1 to find a solution that is optimal or close 
to optimal, but not proven. Proving optimality by using an OPTCR less than 0.01 is 
not recommended since this may result in an excessive solution time and excessive disk 


storage requirements. 
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